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I. 


Introduction 


In order tg better understand and explain the complex, interactive 
effects influencing water movement and water quality in natural systems, 
several mathematical models based on the laws of conservation of mass, 
momentum and energy have been developed (1"7) . Models describing the hydro- 
dynamic and material transport behavior of Mobile Bay have been formulated 
and tested under a National Aeronautics and Space Administration, Environ- 
mental Applications Branch Contract (NAS8-29100) . 

This booklet shows the way in which the model describing water movement 
and tidal elevation is formulated, computed and used to provide basic data 
about the system. Mobile Bay is used in a case study with comments as to 
how the formulations might be expanded or focused to describe other areas 
which qualify under the model restrictions and assumptions. 

The hydrodynamic model, as it will be called throughout the booklet, 
is based on two-dimensional, unsteady flow equations. The water mass is 
considered to be reasonably mixed such that integration (averaging) in the 
d.epth direction is a valid restriction. Convective acceleration, the Coriolis 
force, wind and bottom interactions are included as contributing terms in 
the momentum equations. The equations which makeup the hydrodynamic model 
include the continuity, x-momentum and y-momentum equations (Table 1). 


Table 1, — Mathematical Representation and Operational Modes of 
Physical Models for Mobile Bay 

the 

Nanc 

Eouarion Form 

Results 

Modes 

Concinuicy 

5Qx + 5Qy 
oy 

+ Is « -(R + E) 

at 

Tidal Height 

Tidal Cycle 
Daily Avg 
Monthly Avg 
Seasonal 

:Mon:&ncum 

x**Coaiponifint 

Ms + go M 
be ax 

= KV^cos* - fQQxO'2 
+ Qx(2Wsincj)) 

X- Component of 
Surface Current 

Tidal Cycle 
Daily Avg 
Monthly Avg 
Seasonal 

y- Coaponent 

JSi: + gD 
3t by 

- KV^sint- fQQyD*^ 
+ <3y(2Hsin9) 

y-Component of 
Surface Current 

Tidal cycle 
Daily Avg 
Monthly Avg 
Seasonal 


Results can be calculated for unsteady flow when boundary conditions are 
available as a function of time (dynamic) , or for quasi-steady flow when con- 
ditions are stable for a time period encompassing several tidal cycles or 
longer periods (i.e. weekly, monthly, seasonally or yearly averages). 

The solution of the equations shown in Table 1, applied to Mobile Bay, 
have been used to investigate the influence that river discharge rate, wind 
direction and speed, and tidal condition have on water circulation and holdup 
within the bay. Storm surge conditions, oil spill transport, artificial 
island construction, dredging and areas subject to flooding are other topics 
which could be investigated using the mathematical modeling approach. 
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To understand how the model might be applied to these topics, this 
booklet is subdivided into four parts for the convenience of the reader. 

These are, in order; 

Basic Concepts in Appling the Hydrodynamic Model to a Real System 

• Model Input Requirements 

# A Detailed Illustration; Application to Mobile Bay 
Q Hydrodynamic Model; Program Listing 

Each section will be discussed separately, however, it is advisable that they 
be covered sequentially during the first reading to reinforce the basic con- 
cepts needed to understand and apply the model. 
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II. Basic Concepts in Applying the Hydrodynamic Model to a Real System 

The detailed mathematical development of the model will not be dis- 
cussed here. Reference to "A Hydrodynamic and Salinity Model for Mobile 
Bay," BER Report No. 168-112 (NAS8-29100) , January 1974, by D. 0. Hill and 
G. C, April is suggested for those interested in such a treatment. The model 
equations used throughout this booklet are those mentioned previously 
(Table 1, page 1) , • 

The model system used for illustrative purposes in this section 
resembles that of Mobile Bay. Some features have been added/deleted to 
include features of the model not necessarily applicable to the Mobile Bay 
system. Users of the model should become familiar with the terminology and 
methodology discussed in this section. This information will be particularly 
helpful (a) in specifying the kind and form of data needed to execute the 
model for a specific system (Section III, page 17) and (b) in interpreting 
results calculated by the model (Section IV, page 26). 

A. Real System Data 

Before a mathematical model can be formulated, knowledge of the 
system to which the model is applied must be obtained. This includes both 
formal (recorded) and Informal (experience) data. Some of the sources of 
data are discussed in this section. 

1. Area Maps - Maps showing the geographical and navigational 
features of the system to be modeled are valuable sources of 
data. From these maps (or charts) information about the 
shape and size of the water mass, river entrances, natural 
and man-made passes, marshes adjacent to the water mass, 
sand bars and islands, channels, etc. can be identified and 

:located. 

Tv 

2. Bathemetric Charts - The depth of water in an estuary can be 
found in navigation maps or ocean survey charts. Mean sea 
level is usually chosen as the reference plane on these 
charts. If the source of data uses a different reference 
plane, adjustments to mean sea level must be made. Some 
additional references used in charts are mean high water and 
mean low water. 

3. Tidal Information - Tide heights are taken from appropriate 
tidal charts. Data measured by gauging devices near the 
system boundaries are preferred. However, when the system 
interacts with the open sea, standard tidal periods for the 
open sea can be used. Conversion of the data must be made to 
coincide with the model reference plane - mean sea level. 

The converted data can then be represented by a least squares 
method using a Fourier series of the form; 

HXX = Cl + C2 * COS (C3 * t + C4) 
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where 

HXX is the elevation of water above the reference plane 

Cl is a constant used to adjust the x-axis to a plane 
of symmetry 

C2 is a constant related to the tidal amplitude 

C3 is the tidal frequency 

C4 is the phase angle 

t is the time 

4. Wind Data - Tabulated data and data obtained from wind charts 
are provided by the weather bureau. When the system is large 
(i.e. Mobile Bay) wind conditions at several locations over 
the bay are necessary. This, however, is often impossible 
data to find and good judgement must be used when considering 
wind effects in the model. For periods of relatively steady 
wind direction and speed, no problem is encountered. Like- 
wise, when trend analyses for long periods of time (i.e. 
monthly, seasonally, etc.) are desired, averaging methods can 
be used with good accuracy. However, when wind conditions 
are highly variable over short spans of time, the data fre- 
quency, precision and accuracy becomes critical. 'As a rule of 
thumb, wind speeds below 10 knots are considered to have little 
influence on the hydrodynamics, while winds approaching 15-25 
knots can have pronounced effects on water movement and 
holdup in the system. 

5. River Discharges - River flow information is usually available 
from state geological surveys and/or United States surveys. 

When the system is located in a coastal region where navigation 
is prevalent, data on river flows may also be available in 

U. S. Corps of Engineer files. These data sources may not 
coincide with system boundaries requiring good engineering 
judgement in their interpretation and application. 

B, Model System Features 

Once the real system is known and the area of interest defined, the 
next step in the process deals with the construction of a model system having 
all of the features of the real one. Because the model is governed by com- 
putational laws and rules, the model system will only approximate the real 
system. Care, however, must be exercised to include as many of the important 
features of the real system as is economically feasible. Obviously, the 
real system can be exactly duplicated if the model dimensions are made 
to coincide with real system dimensions. However, computation time would 
become increasingly prohibitive. Once again good judgement must be 
exercised. 
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1. Model Dimensions: Grid Size - The choice of the grid size 

to be used in the model is a function of (a) the real system 
size, (b) the time increment for which a solution is desired 
and (c) the available computer and its costs (Figure 1). 
There is a direct relationship between the incremental time 
for which a solution is desired, model system grid size and 
real system depth dimensions. This is expressed in the form 
of a stability criteria relationship of the form: 


t < 


AS 


v^gH 


max 


where H 

max 


is the maximum water depth encountered in 
the system 


g is the gravitational constant 

AS is the grid dimension 


t is the incremental time for which a solution 

can be specified to insure stability 

The total number of grids should not exceed 2500, Since the 
time increment plays an important role in the total computa- 
tion time, the number of grids selected is often constrained 
to keep computer time reasonable. Thus, final selection of 
the grid size for a system requires justification among the 
following variables : 

a) size of the real system 

b) time increment over which results are desired 

c) economy (computer time), 

2. Water Cell Segments - Water cells on each row are separated 
into segments by land cells encountered on a given row. For 
example, water cells on row 4 in Figure 2 are separated into 
two segments as are the cells on row 9. The water cells for 
most parts of the bay are segmented only once (as in row 7, 
Figure 2.) A water cell segment is labeled by the left most 
(IBNDL) and right most (IBNDR) cell numbers. 

3. Sea Boundaries - Boundaries adjacent to open seas are defined 
as sea boundaries. These boundaries are classified into four 
categories (Figure 3). 

. left sea bounaaries (SL) 

. right sea boundaries (SR) 

. top sea boundaries (ST) 

, bottom sea boundaries (SB) 



Fig. 1 - Grid 



Fig. 2 - Water Cell Segments 
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Fig. 5 - Sea Boundaries 
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Left and bottom sea boundaries are considered to be tidal 
driven boundaries. These are the source of the tidal driving 
force which interacts with the estuary. 

Necessary inputs for sea boundaries include type, total 
number of cells, cell numbers, and tidal conditions which 
drive the model through the tidal cycle. 

4. River Boundaries - River boundaries are boundaries where 

streams enter the system, and/or marsh areas adjoining the 
main water system are located. River discharge cells are 
separated into four categories (Figure 4) : 

^ left river boundaries (RL) 

9 right river boundaries (RR) 

f) top river boundaries (RT) 

^ bottom river boundaries (RB) 

Necessary inputs for river boundaries include type, total 
number of cells, cell number and information relating the dis- 
charge rate into the system. 

River discharge rates are usually modified to reflect 
the influence that each river entrance has on the main water 
movement (i.e, usually expressed in terms of cubic feet per 
second per foot of cell dimension) . An alternate approach is 
available, however, when rivers are located adjacent to marshes. 

River -marsh boundaries are boundaries designed to 
simulate the action occurring in areas which undergo rapid 
changes in stream dimension, velocity and hence discharge rate. 
By including overflow weir equations along the river course 
through a marsh area, drainage into or from the marshes can 
be regulated to produce a known water condition within the 
entrance cell. The marsh area serves as a capacitance element 
which stores river water (Figure 5a, page 11) when tidal con- 
ditions (flood) retard its entrance into the estuary, and 
supplies river water when tidal conditions (ebb) are less 
restrictive. This method also allows a constant and true 
river discharge to be introduced instead of the adjusted 
value required if marsh interaction is omitted. 

River -marsh boundaries are only allowed on the top part 
of the estuary with the weir boundaries parallel to the y-axis. 
Choice of weir heights is a trial and error process depending 
on the tidal influence and size of the system under investi- 
gation. 





Consta 

River 


Discha; 



Fig^ 5a - Rive 









13 


5. Ship Channels - Channels exist in many estuary systems. The 
right side and/or top edge of the water cell nearest the 
channel are used to describe the existence of a channel pass- 
ing through the water mass. The axis (either x or y) of the 
channel, and the depth of the spoil bank are the data 
necessary for the program to describe the ixifluence of the 
channel. For example, channel section A in Figure 6 is 
described by cell 69y and 82y while channel section B is 
described by cell 83x. 

6. Bottom Roughness - Bottom friction also affects flow 
characteristics and requires input of information. Manning's 
equation is used in the program to evaluate bottom frictions. 
Variations of bottom roughness are estimated by altering the 
Manning coefficients which best describes the kind of bottom 
located in each cell. These coefficients usually vary from 
0.010 to 0.050. Marsh areas which have low flow rates may 
have a value close to 0.050. Oyster beds, channels and spoil 
banks may also have relatively high bottom friction 
coefficients. 

7. Passes - The general program has the ability to calculate 
flow rates through vertical and horizontal passes (Figure 7). 
Flow rates through passes which are neither vertically or 
horizontally located can be calculated by arranging imaginary 
passes which have cells in a row or a column close to the real 
pass. A total material balance which involves water through 
the passes and water discharged by rivers over a tidal cycle 
provides a means of verifying the model of a system, 

8. Cell Modifications - Cell modifications are used in the 
calibration phase of the study to adjust calculated results 
in a way that more closely agrees with observed behavior in 
the real system. These modifications are especially needed 

to introduce effects which cannot be represented by the model. 
For example, the channel dimensions are usually small relative 
to the dimensions of the estuary through which it passes. 
Selection of the grid size to affect solution of the model 
equations may seriously neglect a strong local effect caused 
by the channel. When these strong, but unseen (from a model 
point-of-view) effects are observed, cell modifications during 
calibration permits the incorporation of the effects. Two 
ways to introduce modifications include adjustments to bottom 
friction factors, and the inclusion of underwater barriers 
which restrict flow (Figure 8). 

Section II presented some of the terminology and definitions used in the 
establishment of a model system. Section III includes a detailed listing of 
the variables, which, when introduced into the model program, permit one to 
approximate the real system. These variables are given by name, computer 
label and function. Their use will become more clear in Section IV where 
they are used to describe the Mobile Bay system. 
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Fig. 7 - Passes 
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Pig. 0 - Cell Modifications 
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III. Model Input Requirements 

This section lists all the necessary input variables, formats of the 
input variables, aind necessary remarks for defining and using the variables 
to establish a model. The information is presented in the order in which it 
is introduced into the computer. Cross reference of material in Section III 
with the illustration in Section IV is recommended to better understand the 
function of each variable. 

Program Input 

1. Input #1 - Program Control Card 

Format: (1112) Right Justified Integer Number 

Column Variables Remark 

I- 2 IFMT 0, 5 or 6 ; control for width of print field 

0 = default to 6 print positions 

5 = print field width will be 5 print positions 

6 = print field width will be 6 print positions 

3-4 ICORF 0 or 1; control for Coriolis force calculation 

0 = Coriolis force equal to 0 

1 = Coriolis force equal to 0.0007.7 

5-6 ICONV 0 or 1; control convective acceleration 

calculation 

0 = skip convective acceleration calculation 

1 = perform convective acceleration calculation 

7-8 lORIG 0 or 1; control for resistance term calculation 

0 = skip resistance term calculation 

1 = perform resistance term calculation 

9-10 IRMB 0 or 1; does the estuary system have river 

marsh boundaries? 

0 = no river marsh boundary 

1 = perform river marsh boundary calculation 

II- 12 ISPBK 0 or 1; spoil bank calculation for ship channels 

0 = no ship channel 

1 = perform spoil bank calculation for ship 

channels 

13-14 IDSCH 0 or 1; control for discharge rate calculation 

0 = skip discharge rate calculation 

1 = calculate discharge fate at specified 

location 

15-16 IMODF 0 or 1] control for friction factor modification 

0 = use default value of 1 

1 = modify thefriction factors of specified 

cells 


nwEomomjLm of the 

OEIGINMi PA0E IS POOE 
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17-18 

ICK 


0 or 1; read control for velocities from 
previous run 

0 = first run, skip the read statement 

1 = read velocities from previous run 

19-20 

IMNDP 


0 or 1; printing control for manning factors 
and initial depths 

0 = skip printing 

1 = print manning factors and initial depth 

21-22 

INC 


Integer number; print interval 

1,2,3, or ... = print intermediate velocities 
heights, etc. after the specified 
number of increments of time 

Input #2 
Format : 

- Geometrical Information 
(715) Right Justified Integer Number 

Column 

Variables 

Remark 

1-5 

ISB 


Number of total water cell segments; maximum 
number 200. 

6-10 

IROW 


Number of rows, ; maximum number 200. 

11-15 

IV 


Number of cells that require zeroing of 
velocity in y-direction; maximum number 250. 

16-20 

lU 


Number of cells that require zeroing of 
velocity in x-direction; maximum number 250. 

21-25 

IR 


Number of rivers being considered; maximum 
number 20. 

26-30 

NC 


Total number of cells; maximum number 2500. 

31-35 

NCPR 


Number of cells per rows. 

Input #3 
Format : 

- Physical Information 

(6F10.3) Right Justified Real Number with Three Decimal 
Digits or Number with Decimal Point 

Column 

Variables 

Remark 

1-10 

W 


Wind velocity in ft/sec. 

11-20 

THETA 


Wind direction measured from x-axis counter- 
clockwise. 

21-30 

DELS 


Grid size in ft. (square grid). 

31-40 

DELT 


Time increment in sec. 
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41-50 

TIM 

Maximum time for run in hour. 



51-60 

SLADJ 

Sea level adjustment; reference plane is 
mean sea level. 

4. 

Input #4 
Format : 

- Left Boundary Cells of Water Cell Segments 
(1615) Right Justified Integer Number 



Column 

Variables 

Remark 



1-5 

IBNDL(I) 

IBNDL is the cell number of the left most cell 
in a segment. 


6-10 

IBNDL(2) 

Repeat as many times as total number of 

segments . 

5. 

Input #5 
Format : 

- Right Boundary Cells of Water Cell Segments 
(1615) Right Justified Integer Number 



Column 

Variables 

Remark 



1-5 

IBNDR(l) 

IBNDR is the rell number of the right most 
cell in a segment. 


6-10 

IBNDR(2) 

Repeat as many times as total number of 

segments . 

6. 

Input #6 
Format : 

- Cells that Require Zeroing of Calculated Velocity in y- 
(1615) Right Justified Integer Number 

•Direction 


Column 

Variables 

Remark 



1-5 

IBC(l) 

IBC is the cell number of the cell that 
requires zeroing of calculated velocity 
y-direction. 

in 


6-10 

IBC(2) 

Repeat as many times as IV in Input #2. 


7. 

Input #7 
Format : 

- Cells that Require Zeroing of Calculated Velocity in 
x-Direction 

(1615) Right Justified Integer Number 



Column 

Variables 

Remark 



1-5 

IBD(l) 

IBD is the cell number of the cell that 
requires zeroing of calculated velocity 
x-direction. 

in 


6-10 

IBC(2) 

Repeat as many times as lU in Input #2. 
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8. 

Input #8 

- Discharge Rates of Rivers 


Format : 

(8F10.5) 
Digits or 

Right Justified Real Number with Five Decimal 
Number with Decimal Point 


Column 

Variables 

Remark 


1-10 

RB(1) 

Input discharge rate of each river in 
cu.ft./sec. per ft. of grid length. 


11-20 

RB(2) 

Repeat as many times as number of rivers. 

9. 

Input #9 

- Manning 

Coefficients 


Format : 

(26F3.0) 
or Number 

Right Justified Real Number with Zero Decimal Digit 
with Decimal Point 


Colximn 

Variables 

Remark 


1-3 

AMNY(l) 

Input Manning coefficient for each water 
cell segment by segment beginning from first 
segment. 


4-6 

AMNY(2) 

Repeat as many times as number of cells in a 


segment, maximum 26 data per card, continue 
on a new card; start a new card for a new 
segment. 


10. Input #10 - Depths of Water for Water Cells Measured from Reference 
Plane 

Format: (26F3.0) Right Justified Real Number with Zero Decimal Digit 

or Number with Decimal Point 


Column Variables 

1-3 Z(l) 

4-6 Z(2) 


Remark 

Input depth of water of each water cell 
segment by segment beginning from first segment 
If the depth is measured from a reference plane 
other than mean sea level, an adjustment should 
be made by "SLADJ" in Input #3. 

Repeat as many times as number of cells in a 
segment, maximum 26 data per dard, continue 
on a new card; start a new card for a new 
segment. 


If IMODF = 0, skip to Input #14 . 
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11, Input #11 (Optional) - Modification of Friction Factors 
Format : (215) Right Justified Integer Number 


12 . 


13. 


Column 

Variables 

Remark 



1-5 

NMODFX 

Total number of 
the x-dir,ection 

cells which friction factors 
(FX) have to be modified. 

in 

6-10 

NMODFY 

Total number of 
the y-direction 

cells which friction factors 
(lY) have to be modified. 

in 


Input #12 (Optional) - Cell Number and Value of Friction Factor in 

x-Direction for each Cell being Modified 


Format : 

(8(I5,F15,2)) 

- A Right Justified Integer Number and a Right 
Justified Real Number with 2 Decimal Digits 

Column 

Variables 

Remark 

1-5 

KC(1) 

Input the cell number of the cell being 
modified. 

6-10 

VF(1) 

. Input the value of friction factor in 
x-direction for the cell. 

11-15 

KC(2) 

Repeat KC and VF as many times as NMODFX 
in Input #11. 

16-20 

VF(2) 


Input #13 

(Optional) - 

Same as Input #12, Replace FX by FY, x-Direction 
by y-Direction, and NMODFX by NMODFY 


14. 


Input #14 - Starting Time and Total Number of Tide Station 
Format: (F5.2,I5) Right Justified Real Number with 2 Decimal Digits 

and Right Justified Integer Number 

Column Variables Remark 

1-5 TMIL Input the chosen starting time in hours (24 

hour cycls.) for a run. 

6-10 ITSN Input the total number of tide stations in 

the system. 






22 


15. Input #15 - Tide Equation Coefficients and Time Phase 

Format: (5F10.6) Right Justified Real Number with 6 Decimal Digits 

or Number with Decimal Point 


Column Variables 


Remark 


I- 10 Cl(l) 

II- 20 C2(l) 
21-30 C3(l) 
31-40 C4(l) 


Input coefficients of tide equation of the 
form: 

IIXX = Cl + C2 * Cos (C3 * T + C4) for 
Station 1. 


41-50 T(l) 


T is the time variable in the equation. 


I- 10 Cl(2) 

II- 20 C2(2) 
21-30 C3(2) 
31-40 C4(2) 


Start a new card for Station 2 and input the 
coefficients and the time variable as before. 


41-50 T(2) 

Repeat the same for Station 3, 4, ... 
as many times as total number of station. 

16 . Input #16 - Total Number of Boundaries 
Format: (15) Right Justified Integer 


Column Variables 


Remark 


1-5 NTB 


Input total number of river and sea boundaries 
except river marsh boundary. 


Repeat Input #17A, 17B, and 17C as many times as NTB in Input #16 . 


17. Input #17A - Boundary's Nature, Total Number of Cells in the Boundary, 

and Relative Information Required to Perform Calculations 
Format: (A2, 13, FlO.3) Two Alphabetic Characters, Right Justified 

Integer and Right Justified Real Number with Three Decimal 
Digits 


Column Variables 


Remark 


1-2 NATURE (I) 

3-5 NBC (I) 

6-15 RINF(I) 


Boundary natures are classified into eight 
categories (See Fig. 3 and 4). Each category 
requires its relative information to perform 
the calculation. 

1. SL - Left sea boundary 

2. SB - Bottom sea boundary 

3. SR - Right sea boundary 

4. ST - Top sea boundary 
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The relative information for a sea boundary 
is the related tide station number, 

5. RL - Left river boundary 

6. RB - Bottom river boundary 

7 . RR - Right river boundary 

8. RT - Top river boundary 


The relative information for a river boundary 
is the river discharge rate in cu.ft./sec. per 
feet of cell dimension. 

NBC is the total number of cells in the boundary 


Input #17B - Cell Numbers of the Boundary 
Format: (1615) Right Justified Integer 


Column 

1-5 

6-10 


Variables 

INBC(I,1) 

INBC(I,2) 


Remark 

Input the cell numbers of the boundary 
beginning from left to right or bottom to top. 

Repeat as many times as NBC (I) in Input #17 A. 


Input #17C (Optional) - Depth of Water Adjacent to Left Sea Boundaries 

or Bottom Sea Boundaries 

Format: (16F5.0) Right Justified Real Number with Zero Decimal Digit 

or Number with Decimal Point 


Column Variables 

1-5 ZB(1) 


6-10 


ZB(2) 


Remark 

Input the depth of water left to the left sea 
boundaries or below the bottom sea boundaries, 
cell by cell. Skip this part if the boundary 
type is not SL or SB . 

Repeat as many times as NBC(I) in Input #17 A. 


18 . 


If IRMB = 0, skip to Input #21 , 

Input 7#18 (Optional) - River Marsh Boundary Bottom Cells 
Format: (215) Right Justified Integer Number 


Column Variables 


1-5 

6-10 


IRMBB(l) 

IRMBB(2) 


Remark 

Input the two cell numbers of the river marsh 
boundary bottom cells. 
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19. Input #19 (Optional) - River Marsh Boundary Top Cells 

Foirmat: (215) Right Justified Integer Number 

Column Variables Remark 

1-5 IRMBT(l) Input the .two cell numbers of the river 

6-10 IRMBT(2) marsh boundary top cells. 

20. Input #20 (Optional) - Depth of Weir Along River -Marsh Boundary 

Format; (FlO.5) Right Justified Real Number with 5 Decimal Digits or 

Number with Decimal Point 

Column Variables Remark 

1-10 ZRB Input the depth of weir in feet along the 

river marsh boundary. The river marsh boundary 
should locate at the top part of the system 
with the boundary weir parallel to y-axis. 


If ISPBK = 0, skip to Input #23 . 


21 . 


22 . 


Input #21 (Optional) - Total Number of Cells of Spoil Banks Along the 

Ship Channel 


Format ; 

(15) - Right Justified Integer Number 

Column 

Variables 

Remark 

1-5 

NSBC 

Input the total number of cells of spoil banks 
along the ship channel. 

Input #22 
Format ; 

(Optional) - Cell Number, Axis of the Spoil Bank and Depth of Spoil 
Bank for Each Cell Along the Ship Channel 
(8 (15 , Al, F4. 1) ) Right Justified Integer Number, Alphabetic 

Character, and Right Justified Number with 1 Decimal Digit 

Column 

Variables 

Remark 

1-5 

INSBC(l) 

Input the cell number. 

6 

lAXIS(l) 

Input the spoil bank axis either x or y. 

7-10 

ZCH(l) 

Input the depth of spoil bank in feet for 
the cell along the ship channel. 

11-15 

16 

INSBC(2) 

IAXIS(2) 

Repeat as many times as NSBC in Input #21; 
maximum eight data group a card. 

17-20 

ZCH(2) 




If IDSCH = 0, skip Input #25 . 


23. Input #23 - Total Number of Passes Where Discharge Rates will be 
Calculated 

Format: (15) Right Justified Integer 


Column 

Variables 

Remark 

1-5 

LDSCH 

Input the total number of passes where dis 
charge rates will be calculated. All the 
cells at a location should be either on a 
row or a column. 

Input #24 
Foinnat : 

- Total Number 
(16(14, Al)) 

of Cells at a Pass, and the Flow Direction 
Remark 

Column 

Variables 

1-4 

NCDSCH(l) 

Input total number of cells at the pass. 

5 

lAXIS(l) 

Input the flow direction either x or y. 

6-9 

10 

NCDSCH(2) 

IAXIS(2) 

Repeat as many times as LDSCH in Input #23 
maximum 16 data group a card. 

Input #25 
Format : 

- Cell Numbers at Each Pass 

(1615) Right Justified Integer Number 

Column 

Variables 

Remark 

I- 5 
5-10 

II- 15 

ICDSCH(1,1) 

ICDSCH(1,2) 

ICDSCH(1,3) 

Input the cell numbers at Pass 1 from left 
to right or from bottom to top; repeat as 
many times as the number of cells in the 
pass . 

1-5 

5-10 

ICDSCH(2,1) 

ICDSCH(2,1) 

Start a new card for Pass-2 and input the 
cell numbers as before. 


Repeat the same for Pass 3, 4, ... 

As many times as total number of stations. 
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IV. A Detailed Illustration: Application to Mobile Bay 

This section illustrates the use of the general program beginning from 
data collection to model verification as applied to the Mobile Bay system. 

In order to demonstrate all of the model features, the real system behavior 
has been modified to include elements which do not necessarily exist in the 
bay at this time (i.e. the existence of an island near the Theodore area). 
Figures are inserted at appropriate places to help explain the procedure. 

The format of this section is established to coincide with those in Section 
II and Section III. 

A. Mobile Bay data: 

1. Map and Charts; / 

. Obtain a navigational chart of Mobile Bay and vicinity area 
(Figure 10) published by the U. S. Coast and Geodetic Survey 
(^1^1266 ) 3 location map of South Alabama (Figure 9) from 

the Geological Survey of Alabama. 

Locate the following system features; rivers, islands, 
channels, passes, sand bars, head lands, river marsh area, etc. 

Sketch the shape of the water system selected for modeling 
and measure the horizontal and vertical distances which enclose 
the selected system. 

2. Bathemetry; 

Bottom elevations (Z(I)) from mean water are taken from the 
navigational chart (Figure 10) . Adjust the values to the 
reference plane - mean sea level. The difference between mean 
sea level and mean water of the chart is 1.8 (SLADJ) supplied 
by the Corps of Engineers, Mobile District. 

3. Tidal Information; 

Tide charts (Figure 11 to Figure 15) were provided by the 
Corps of Engineers, Mobile, Alabama for Mobile State Docks, 
Dauphin Island-Gulf, Cedar Point, Bon Secour, and Great Point 
Clear. 

Tabulate the time and tide height data for Dauphin Island-Gulf, 
and Cedar Point from the charts. Fit the data by least square 
to the equation of the form (page 30-31): 

HDI = Cl + C2*Cos (C3 * t + C4) 

where HDI = tide height 
t = time 

Cl, C2, C3, C4 are coefficients 
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The following are the least-square fit equations: 

Dauphin Island - 

HDB(I) = 1.090 + 1.295 * cos (0.004188**0.0567.4) 

Cedar Point - 

HDB(I) = 1.089 + 1.177 * cos (0.004188*-0. 0032453) 

4. Wind Data: 

Recorder charts (Figure 16) with tabulated data were provided 
by the National Oceanic and Atmospheric Administration- 
Environmental Data Service at Ashville, North Carolina. 

Calculate the average wind velocity in feet per second 
and direction measured from x-axis counter-clockwise from 
the charts, and use these as input data. 

W = 5.60 knots 

= 9.46 feet/second 

THETA = 70° 

5. River Flows: 


Flow information (Figure 17 and Figure 18) for the Alabama 
River at Claiborne and the Tombigbee River at Coffeeville 
was provided by the U. S. Geological Survey, Montgomery, 
Alabama. Calculate the average discharge rate in cubic 
feet per second for each river and convert these values to 
flow rates per grid unit length. The flow rate of the 
Alabama River at Claiborne, Alabama from the chart is: 


44000 ft3/ sec 
" 6561.68 ft 


6.706 ft^/sec 


The flow rates for the other rivers are considered zero 


this example. 


td 

It 

O 

(Dog River) 

RB(3) = 0 

(Fowl River) 

RB(4) = 0 

(Bon Secour River) 

RB(5) = 0 

(Fish River-Weeks Bay) 

RB(6) = 0 

(Tensas River) 
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Fig. 16 - Typical Weather Data for the Mobile, Alabama Area. 
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B. Model System Features of Mobile Bay 
1. Grid Size 

Superimpose a two-dimensional grid over the bay area (see 
Figure 19). Label the cells and differentiate water and land 
cells. Grid spacings in x and y directions are equal and 
chosen as tw'j/ kilometers (6561.68 feet (DELS)). There are 
21 cells (NCPR) per row and 38 rows (IROW) , giving 798 
cells (NC) . Since the deepest water in the bay is 40 feet, 
this sets the upper limit of the incremental time for cal- 
culations at 

_ . AS 

^ /'^'g'Hma7 

At < 130 sec. 

DELT = 120 i^ chosen as program input. 

After setting the grids, record the depth of water from mean 
water in each water cell from the navigation chart. 

Cell No. 


1 
1 

2 

3 

4 

22 

23 


97 

98 

99 


787 

788 

789 


Depth 

Z(I) 

40 

40 

30 

39 


37 

34 


3 

6 

5 


1 

1 

10 
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2. Water Cell Segments (Figure 20) : 

Label each water cell segment by the left most boundary cell 


(IBNDL) and the 

right most boundary 

cell (IBNDR) 

• 

Row No . 

L. Cell No. R. Cell No. 


N 

IBNDL (N) 

IRNDR(N) 


1 

1 

11 


2 

22 

32 


3 

43 

53 


4 

64 

70 


5 

89 

94 


6 

97 

99 


7 

109 

124 


Sea Boundaries 

(Figure 21): 



Label the sea boundaries by SL, SB, 

SR, or ST. 

Locate tide 

stations, set up 

the necessary data. 



S. Boundary No. 

Type No. of Cells Rel. 

Inf n . 

I 

NATURE (I) NBC (I) RINF(I) 

1 

SB 10 


1. 

2 

SL 1 


2. 

3 

SL 4 


1. 

4 

SR 3 


1. 

sequence not 




important 




S . Boundary No . 

No. of Cells 

Cell Nos . 


I 

J = NBC(I) 

INBC(I) 


1 

10 

2 




3 


2 

1 

149 


4 

3 

11 




32 




53 


For SL and SB, 

• • 

also supply boundary depth, i.e. 

the depth of 

cells left to or below the boundary 
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Pig. 20 - Water Cell Segments 
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Boundary No, 


No. of Cells 


Boundary Cell Depth 


I 

1 


J = NBC (I) 
10 


ZB (I) 
40. 
30. 
39. 


1 

4 


8 . 

42. 

40. 


4. River Boundaries (Figure 22); 

Locate the river discharge cells, and the river-marsh area. 
Separate the boundaries into their ow categories. For river 
boundaries, tabulate the number of cells and flow rates. 

Small streams may be neglected. 


R , Boundary No , 

1 
1 

2 


'h2± 

NATURE (I) 
RT 
RL 


No. of Cells 

NBC (I) 

1 

1 


Flow Rate 

RB(I) or U(I) 
6.706 
Q.iez 


For river -marsh boundaries, locate the river -marsh bottom cells 
and top cells. 

IRNiBB(l) = 595, IRMB(2) - 599, IRMBT(l) = 784, IRMBT(2) = 788. 

5. Ship Channel (Figure 23): 

Locate the ship channel and label the cells along the channel. 
Depths of spoil banks were obtained from the Corps of Engineers, 
Mobile, Alabama, Tabulate the cell numbers, channel axes, and 
depth of spoil banks of each cell along the channel. 


Total No. of Cells 


NSBG 

22 

Cell No. 

Channel Axis 

Depth of Spoil Sauk (Ft) 

INSBC(I) 

lAXlS (I) 

2CH(I) 

153 

y 

9. 

174 

y 

7 . . 

195 

y 

7. 

196 

X 

9. 






























































Fie. 25 - Ship Channel 





























































































Bottom Roughness (Figure 24) ; 


Manning coefficients which are used to calculate the bottom 
friction vary from 0.010 to 0.050. A coefficient of 0.05 is 
used in the marsh area to simulate the low flow rates 
expected to this area. Values within the bay range from 0.010 
to 0.018. The coefficients of water cells are read in segment 
by segment. 

3 

Segment Manning Coefficients x 10 


Z(I) 


1 

15. 

15. 

15. 

15. 

15. 

15. 

15. 

15. 

15. 

15. 

2 

• 

15. 

15. 

15. 

15. 

15. 

15. 

15, 

15. 

15. 

15. 

• 

• 

78 

15. 

50. 

50. 

50. 

50. 

15. 






Passes (Figure 25): 

Babel the passes by number. Count the total number of cells 
along each pass. Decide the flow direction either x or y 
through each pass. 

Total No. of Passe s 
LDSCH 
2 

Flow Direction Cell Nos. (Each Pass) 
lAXIS(K) ICDSCH(K,J) 

y 89 

90 

91 

X 149 


Cell Modifications: 

a. Zeroing of velocities: 

Cells with land boundaries or sand bars on top edges 
have zero flow rates across the top edges. The same 
implication applies to cells with land boundaries on the 
right side. Count total number of cells that require 
zeroing of Vp and Up (Figure 26 and Figure 27) 
respectively, and record the cell numbers. 


Total No. Cells (Each Pass) 
NDSCH(K) 

2 3 



?lg. 24 - Manning Coefficionta 
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Fig. 25 - PasBOD 
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Total No. of Cells, Vp = 0 Total No. of Cells, Up = 0 


Cell No., Vp=0 
IBC(I) 

25 26 64 65 66 67 etc. 

Cell No.. Up=0 
IBC(I) 

70 94 298 -320 99 124 


b. Modifying Friction Factor; 

Pilot runs are necessary to obtain friction factor 
modification information. The following modification 
is used to restrict flow in the upper part of the bay,. 

Count total number of cells which FXs have to be 
modified and which FYs have to be modified, respectively. 
Record cell nos. (Figure 28), 

Total No. Cells , FX modified Total No. Cells , FY modified 
NMODFX NMODFY 

8 8 ' 

Cell No. & Value of FX KC(I ), VF(I) 

362 1.1 384 1.1 '406''' 1.1 


Cell No. Sc Value of FY KC(I) , VF (I) 

383 1.1 .405 l.T'”' 427" 1.T"™. .. 

C. Model Input for Mobile Bay (Refer to page 17 for variable 
description and pages 26-51 for Mobile Bay values.) 

1. Input #1 “ Program Control Card 

IFMT, ICORF, ICONV, lORiG, IRL-ffi, ISPBK, IDSCH, IMODF, ICK, 
IMNDP, INC. 


2. Input #2 - Geometrical Information 

12.34bb7a9l)i234667b90 5, 23‘|507M9U1£345670901 2345 2^4567690 12 

0 1 1 1 i 1 i 1 0 0 4 

l23'4bfa7B9U12345fa7tig0l2345G7a90i2v34 567B901234567t)99l23‘ 567090 1234 56709U12 
39 30 30 30 6 790 21 



! V TJTg 
hi io r'^OR 
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Input #3 Physical Information 
W, THETA, DELS, BELT, TIM, SLADJ 


53 


li:34b67a90ia3H5670901i;445678901234567B90l234567B90l23i 567090123456789012 
0. 0. 6561.68 120. 25. 0. 0.0 


4. Input #4 “ Left Boundary Cells of Each Segment 
IBNDL 


l.i3456789Ul2345fa7og0l234567H9ui234567R901234567090l2345678y0 12345678901234567890 


2 

23 

44 

65 

89 

97 

109 

129 

149 

170 

l'J2 

213 

235 

2Sb 

277 

29fl 

319 

340 

3oi 

302 

403 

425 

446 

467 

489 

510 

5.)1 

552 

574 

595 

616 

637 

658 

b79 

700 

721 

742 

763 

704 











5. Input #5 - Right Boundary Cells of Each Segment 
IBNDR 

l234bb789U12345670901234567fl9ui234567090l2345b7090l2345670y0123456789Ul234b67HyO 


11 

32 

53 

70 

94 

99 

124 

145 

166 

186 

2U6 

226 

247 

26b 

20 b 

306 

327 

663 

340 

684 

3bd 

705 

389 

726 

411 

747 

432 

768 

454 

709 

474 

495 

516 

5.S7 

558 

579 

6U0 

621 

642 


6. Input #6 - Cells that Require Zeroing of Calculated Vp 
IBC 

I234567a9012345b7090i2345b70901234567a90l2345670901234567890i23456?89ul234b67ug0 


25 

2b 

64 

65 

66 

67 

67 

50 

51 

52 

!)3 

54 

5 b 

5b 

57 

50 

59 

247 

60 
26 b 

bi 

28b 

62 

348 

149 

454 

2lo 

246 

403 

467 

552 

575 

576 

577 

578 

16b 

10b 

206 


54 

7. Input #7 - Cells that Require Zeroing of Calculated Up 
IBD 

1^34&678901234b676g01234567B9Ui23456709012345fa7B90123'>5678901234567090l234b67b90 


70 

94 

99 

124 

145 

16b 

186 

206 

226 

247 

2o6 

286 

298 

306 

320 

327 

340 

360 

3b9 

411 

432 

454 

474 

495 

516 

537 

5'.)Q 

b/9 

600 

621 

642 

663 

604 

705 

726 

747 

768 

789 












8. Input #8 - Discharge Rate of -Rivers 
RB 

1234bt>7a901234567a901234567n90123456.78901234567890l2345fa78901234567fl9Ul2 
6.70b60 U.O 0.0 0.0 0.0 0.0 

9. Input #9 - Manning Coefficient 


I<i34ii67890 1234567690123456789012345678901234567890 1234 567090 123456789012 


15 

15 

15 

15 

15 

15 

15 

15 

15 

15 








15 

15 

15 

15 

15 

lb 

15 

15 

15 

15 








15 

15 

15 

15 

15 

15 

15 

15 

15 

15 








15 

15 

15 

15 

15 

15 












18 

15 

15 

15 

18 

18 












10 

18 

18 















10 

16 

18 

18 

18 

15 

15 

18 

10 

16 

10 

10 

10 

18 

10 

18 


18 

15 

15 

18 

15 

18 

18 

18 

15 

15 

15 

15 

15 

15 

15 

15 

18 

10 

18 

18 

18 

18 

18 

15 

10 

18 

lb 

10 

10 

18 

18 

iO 

16 

18 

18 

18 

18 

18 

18 

15 

16 

16 

18 

lb 

10 

18 

16 

18 

18 

16 

18 

i8 

IB 

18 

18 

15 

18 

16 

18 

18 

Id 

10 

10 

18 

18 

18 



18 

16 

18 

18 

15 

18 

16 

18 

16 

16 

18 

18 

16 

18 




18 

18 

18 

15 

16 

18 

10 

18 

10 

18 

18 

10 

18 





18 

18 

18 

15 

18 

18 

18 

10 

18 

la 

18 







18 

18 

16 

18 

15 

18 

10 

18 

10 

10 








18 

18 

18 

18 

15 

18 

18 

10 

18 









18 

18 

18 

18 

lb 

18 

18 

16 

10 









18 

l8 

18 

18 

15 

16 

10 

18 

16 









18 

18 

10 

18 

15 

18 

18 

18 










18 

18 

18 

10 

lb 

18 

18 

18 










10 

lo 

18 

15 

16 

18 

18 

18 

18 









18 

18 

15 

18 

lU 

18 

18 

16 










lO 

10 

10 

10 

18 

18 

18 

16 

10 









10 

10 

10 

10 

18 

Id 

18 

23 










12 

15 

10 

18 

18 

18 

18 











12 

15 

18 

10 

Id 

18 

18 











12 

15 

18 

18 

10 

18 

18 











12 

lb 

18 

18 

18 

10 

18 











15 

18 

18 

18 

18 

10 












15 

50 

50 

50 

50 

15 












15 

50 

50 

50 

50 

15 












15 

50 

50 

50 

50 

15 












15 

50 

50 

50 

aO 

15 












15 

50 

■50 

50 

50 

15 












15 

50 

50 

50 

50 

15 












15 

50 

50 

50 

50 

lb 












15 

50 

50 

50 

50 

15 












15 

50 

50 

50 

50 

15 












15 

50 

50 

50 

50 

15 
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10 . Input #10 - Depth of Water Measured from Reference Plane 

l«;3^b67a9ul23^D67b90l2A^5fa7fl9ui234567890l234567ti90^23.♦567090l2J456789•0^2 
<♦0 30 39 21 21 30 35 37 40 37 

34 39 2B 17 30 30 27 30 39 40 

35 21 9 9 *;U 30 9 30 35 40 

29 6 11 9 23 30 

5 b 30 20 ib 5 

3 6b 

4 5 9 14 17 16 12 3 7 7 9 7 a 4 5 4 

4 0 10 10 13 12 12 12 12 12 10 9 7 0 8 0 5 

10 10 11 12 12 12 11 10 10 10 11 11 10 7 9 8 8 j 

8 a 10 12 11 11 10 10 10 10 11 11 11 7 9 8 4 

7 10 11 iO 11 10 10 10 10 lo 12 12 12 a 3 

3 10 11 11 12 11 11 10 11 10 10 9 7 3 

7 10 10 11 10 11 10 10 10 9 85 1 

7 9 10 12 iO 11 11 10 10 8 2 

5 9 10 11 10 11 11 10 7 

6 10 10 9 12 11 11 10 b 

8 10 10 9 11 11 10 10 4 

4 10 10 9 11 12 9 7 4 

6 10 10 10 11 10 94 

8 9 9 10 il 10 11 6 

3 7 811111111 9 5 

7 a 9 10 10 11 10 9 

‘ 7 8 9 10 10 10 9 10 5 

666667 0 8 

6 10 9 10 9 a 6 

2 10 10 10 9 7 5 

2 8 10 9 8 b 10 

2 8 3 9 6 4 3 

0 2 4 3 1 6 

8 1 1 1 1 fa 

10 1 1 1 1 10 

10 1 I 1 1 10 

10 1 1 1 1 10 

10 1 1 1 1 lU 

10 1 1 1 1 10 

10 1 1 1 1 10 

10 1 1 1 1 10 

10 1 1 1 1 10 

10 1 1 1 1 10 


11. Input #11 - Modification of Friction Factor 
NMODFX, NMODFY 


Input #12 -• Cell Number and Values of Friction Factors for FX 
KC, OF 


Input # 13 - Cell Number and Values of Friction Factors for FY 
KC, OF 


1234bb789ol234 5fa789Ql234Sb7R90123456709012345fa789012:i456789012345678901234567ti90 

6 b 

362 1.1 364 1,1 406 l-i 428 1.1 450 1«1 4/2 1.1 494 1.1 bib 1.1 

383 1.1 40b 1.1 427 1.1 449 1.1 471 1.1 493 1.1 bl5 1.1 537 1.1 




OP TiiS 



56 


12. Input #14 - Starting Time and Number of Tide Stations 
TMIL, ITSN 

Input #15 - Tide Equation Coeffieients and Time Phase 
CA, CB, CC, CD, TD 

l*:3‘46b7a9oi2’3456789012445b7H9ul234567890la3<t567a90l2:i't567U901EiM!j67fl9Uia34567890 
13. b 3 

1.03042 1.24243 0.004188 0.0867114 210. 

1. 0891b 1.1/665 0.00418R -0.003245 120. 

1.1B329 1.306028 0.004188 -0.182190 60. 


13. Input V/^16 - Total Number of Boundaries 
NTB 


12345678901234 567890123456789ui2345678901234567890l2:54567890123456789U12 
7 


14. Input #17 - Boundary's Nature, Total Number of Cells in the 
Bbundary, and Relative Information 

NATURE, NBC, RZNF 

Input #19 - Depth of Water Adjacent to Left or Bottom Sea 
Boundaries 
ZB 


I234bb789ul2345b709012345b789012345678901234567890i2.54567890123456709012 

S8 10 1. 

2 3 4 5 6 7 8 9 10 11 

40. 30. 39. 21. 21. 26. 35. 37. 40. 37. 

SL 1 2. 

149 

8.0 

SL 4 1. 

2 23 44 65 

42. 40. 3b. 36. 

SH 3 1 . 

11 32 63 

RL 1 .7b2 

467 

RT 1 6. 70 bo 

784 
RT 1 

789 


0 . 
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15. Input #20 - River Marsh Boundary Bottom Cells 
IRMBB 

Input #21 - River Marsh Boundary Top Cells 
IRMBT 

Input #22 - Depth of Weir Along River Marsh Boundary 
ZRB 

i^34bb7B9Uia345676901^34567H90i234567890123456709(J1234567890l2i456789Ui^ 

59b 599 

784 788 

1*8 

16. Input #23 •= Total Number of Cells of Spoil Banks Along Ship 

Channel 

NSBC 

Input #24 - Cell Number, Axis of the Spoil Bank in Each Cell, 

Depth of Spoil Bank 

INSBC, I AXIS, ZCH 

i234bb789ul234b6709bl2345b7890123456789012345b7a9012345678901234b6709Uia34b67b90 


153Y 

9. 

1/4Y 

7. 

A95Y 

7. 

196X 

9. 

217Y 

9. 

?.o8Y 

8 • 

259Y 

9. 

280Y 

31UY 

b. 

322Y 

5. 

343Y 

5. 

364Y 

6. 

385Y 

6 • 

4u6Y 

V. 

427Y 

5. 

448Y 

469Y 

3. 

49UY 

0. 

bllY 

0. 

532Y 

0. 

553Y 

0. 

574Y 

0. 





17. Input #25 -■ Total Number of Passes Where Discharge Rates will 
be Calculated 

LDSCH 

Input #26 - Total Number of Cells at Each Pass, lAXlS of Flow 
Direction 

NCDSCH, TAXIS 

Input #27 - Cell Number S at Each Pass 
ICDSCH 

I234b678901234b67890 1234567890123456789012345676901234567890123456/89812 
2 

3Y IX 
89 90 91 

149 



Tr 300J)a OELT; 2.00 MIL. T= 18.50 


not heisht *i stltion i i » -= 
TIDE height AT STATION ( 21 x 
TIDE HEIGHT AT STATION « 3» r 
VOLUME FLOW RATE AT LOCATION I 
VOLUME FLOW RATE AT LOCATION f 


Tr 360.00 DELIr 2.00 

TIDE HEIGHT AT STATION I II r 
TIDE HEIGHT AT STATION L2I - 
TIDE HEIGHT AT STATION I 31 = 
VOLUME FLOW RATE AT LOCATION I 
VOLUME FLOW RATE AT LOCATION I 


T= 420'.00 DELT- 2.00 

TIDE HEIGHT AT STATION (H r 
TIDE HEIGHT AT STATION 1 21 r 
TIDE HEIGHT AT STATION ( 31 = 
VOLUME FLOW SATE AT LOCATION » 
VOLUME FLOW PATE AT LOCATION i 


Tr 480.00 OFLTr 2.00 

TIDE HEIGHT AT STATION ( II r 

TIDE height at station I 21= 
TIDE HEIGHT AT STATION e 31 = 


VOLUME 

VOLUME 

FLOW 

FLOW 

rate 

RATE 

AT LOCATION ( 

AT location t 

RIVER - 

1 

FLOW 

RATE r -44UUU. 

RIVER - 

2 

FLOW 

RATE = 

R1VF,R - 

3 

FLOW 

RATE r 

river 

4 

FLOW 

RATE = 

river - 

S 

FLOW 

RATE r 

RIVER - 

6 

FLOW 

rate r 


.30 

.87 

1.50 

II = -2379:iQ.91 CU.FT./SEC 

21 = 138607.38 CU.Ft./SEC 


MIL. Tr 19.50 

i.U8 

.59 

1.17 

II r -532609.66 CU.FT./SEC 

21 r 86840.43 CU.FT./SEC 


MIL. Tr 20.50 

-.09 

.34 

.85 

Ur -6S6987.59 CU.FT./SEC 

21 = 28599..07 CU.FT./SEC 


MIL. T= 21. 5U 

-.19 

.14 

.54 

II = -678752.15 CU.FT./SEC 

21 r -31559.35 CU.FT./SEC 


00 CU.FT./SEC 


DO CU.FT./SEC 


00 CU.FT./SEC 


DO CU.FT./SEC 


00 CU.FT./SEC 


DO CU.FT./SEC 


a 


Program Output: 





HCI6HT FBOH BEFERCKCE PLANE TO SURFACE 


37 

36 

35 

31 
33 

32 
31 
30 
29 
2 ? 
27 
26 
25 
2 * 
23 
22 
21 
20 
19 
13 
12 
16 
15 
!<• 
13 
12 
11 
10 

9 

a 

7 

6 

5 

3 

Z 

1 


1.07 
1.09 
1.11 

1.08 
1.07 

■ 1.08 
1 . 07 
1.04 
1.02 
. 97 . 33 

1 . 00 1 . 00 

l.ro .98 .93 

1 . 20 1 . 01 , 37 

.97 .95 

■' .93 

.59 . 59 .E 4 

.57 .60 .58 

.59 .61 .62 

.60 .57 .53 


1.37 
1.30 
1 . 17 
1.05 
.95 
.89 
.89 
.91 
.93 
.96 
.97 



.94 

.97 


.93 

.97 


• 97 

.98 


.99 

.99 

.97 

• 98 

1.00 

1.00 

1.01 

1.03 

1 . 03 

1.06 

1.06 

1 . 03 

1.09 

1.13 

1.10 

1 . 13 

1.14 

1 . 12 

1 . 13 

1 . 14 

1.10 

1.12 

1 . 12 

I. 08 

1 . 10 

1 . 11 

1 . 08 

1 . 09 

1.09 

1 . 07 

1.07 

1.06 

1 .D 5 

1.05 

l.US 

1 . 03 

1.03 

1.03 

1.02 

1 .U 3 

1 . 03 

1 . 01 

l.Ot 

1 . 03 

.99 

.99 

1 . 01 

.36 

• 98 

.99 

. 9 <» 

.96 

.99 

.88 

.90 

.91 

. 76 

.80 

.84 

.61 

.68 

.72 

.63 

.64 

.63 

. 6 « 

.66 

.61 

.55 

.59 

.55 


.00 

.00 

.on 

.00 

.00 

.00 

.qu 

.00 

.00 

. 01 ) 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

. 00 

.00 

.00 

.00 

.00 

.00 

.97 

. 98 

.98 

.96 

.97 

.97 

.97 

.98 

.98 

.98 

.98 

.99 

.98 

.39 

1.00 

1.01 

1.03 

1.08 

1.06 

1 . 12 

1 . 14 

1 . 14 

1.14 

1 . 16 

1.15 

1.15 

1.16 

1.16 

1.16 

1 . IS 

1.17 

1.17 

I. IS 

1.17 

1.17 

1 . 16 

1.15 

1.15 

1 . 15 

1.12 

1.14 

1.13 

1.09 

1 . 11 

1.14 

1.05 

l.OS 

1 . 11 

1.02 

1.06 

1.08 

1.02 

1.05 

1.07 

1.03 

1.02 

1.06 

1.02 

1.03 

1 . 04 

1.02 

1.06 

1.02 

.99 

1.04 

1 . 06 

1.01 

.97 

1.03 

.92 

.97 

.98 

.64 

.53 

.56 

.62 

.59 

.59 

. 57 

.57 

.58 


.00 

.59 


.00 

.55 


.00 

.51 


.00 

.50 


.00 

.50 


.00 

.51 


.00 

. S 3 


.00 

.55 


.00 

.59 


.00 

.65 


.91 

.72 


.96 

.95 


.98 

.98 


1.00 

1.04 


1.06 

1.08 


1.10 

1.11 


1.14 

1.16 

1.19 

1.16 

1.18 


1.17 

1.20 


1.16 



1.15 



1.14 

1 . IS 


1.13 

1.13 


1.14 

1 . 13 


1.13 

1.13 

1 . 13 

1.11 

1.12 

1 . 11 

1.09 

1 . 09 

1 . 12 

1.07 

1 . 07 

1 . 09 

1.06 

1.07 

1. 07 

1.08 

1.07 

1 . 03 

l.UG 

1.09 

1.07 

l.US 

1.07 

1.07 

1.05 

1.09 

i.n 


1 . 16 

.59 

.59 

.56 


1.10 


1 . 10 

1 . 03 

1.05 

1. 11 

1 . 09 

1.09 

1.11 

1 . 11 

1.13 

1.09 

1 . 11 

1 . 16 

1.09 

1 . 11 

1 . 14 

I.'IO 

1 . 10 

1.08 

1 . 13 

1 . 12 

1.05 

1 . 19 

1.20 



1.17 



1.17 

1.17 


1.14 

1.13 

1 . 10 

1^05 

1.03 

1.02 

.97 

. 96 

.97 


Ui 

VO 



velocities IN X OIOECIION 
39 






.U 

.0 

.U 

.0 

>u 

.0 









.0 

.U 

..U 

'.0 

»u 

>.o 









.n 

.0 

.0 

.0 


>0 









.-0 

.0 


r.O 

.0 

.0 









.0 

.0 

^0 


.0 

.0 










.□ 

.0 

’.0 

«0 

'.0 









,.o 

■.0 

^0 

'.0 

-o 

so 









.0 

i.0 

.U 

>.0 


to 









.0 

.0 

.0 

..a 

.0 

‘.0 









*.n 

.0 

„U 

J.0 

,u 

-.0 









to 

.2 

.2 

.9 

1.6 

3.2 








;.o 

.8 

.3 

.9 

1 .5 

2.6 

2.3 








.0 

.2 

.2 

UU 

1.6 

1.9 

2.2 








'.0 

.5 

.2 

1.1 

.4 

-.2 

- 1.1 








-0 

.7 

.2 

.5 

-.8 

- 2.1 

-.1 







.8 

1.7 

1.3 

.1 

-.4 

- 1.6 

.5 

.4 







;.u 

l.i 

.9 

-.1 

- z.a 

.9 

l.S 

1.4 

-.2 






i.O 

.6 

.S 

-.5 

2.4 

2.7 

3.1 

2.6 






.0 

.9 

.3 

- 1,7 

-.2 

2.1 

2.6 

3.4 

2.6 






.0 

.5 

- 1.2 

.4 

-.2 

.9 

1.4 

1.4 







.0 

-.8 

.8 

.3 

-.3 

-.6 

-.1 

-UO 







.n 

-.1 

.1 

-.1 

-.4 

- 2^0 

- 3.0 

- 3.9 

- 2.1 






•.0 

-.1 

.0 

-.4 

-.3 

- 1.1 

- 2.2 

- 2.3 

- 1.3 






.a 

;.D 

.4 

.4 

2.0 

.2 

- 1.0 

- 1.7 

- 1 . 3 ' 






-0 

.9 

.8 

.7 

2.8 

.5 

-.7 

- 2.4 

- 2.9 

-1 .4 





-n 

.8 

.8 

.3 

-.1 

- 1.6 

- 2.1 

- 3.1 

- 3.7 

- 3 1.0 

-.9 




• n 

.8 

.8 

.6 

.0 

-.3 

- 1.1 

- 1.8 

- 2.3 

- 2.5 

- 2.0 

-1 .3 

-.3 


i.n 

.8 

.5 

.2 

.1 

-.9 

- 1.2 

- 1.3 

- 2,.0 

-1 .7 

- 1.3 

- .7 

- .1 

.0 

T 

.3 

.0 

‘.0 

-.0 

-.3 

-.3 

-.6 

- 1.3 

- 1.2 

- l.l 

-1 .5 

- .5 

-1 .9 

2.9 

3.0 

1.5 

-.1 

-.6 

.6 

-.0 

.3 

-.3 

- .4 

-.4 

- .8 

-1 .0 

2G.2 

IB. 9 

12.1 

S.9 

-!.0 

1.1 

1,3 

.8 

.4 

,4 

.5 

- .4 

-I .0 

- .9 

.0 

2.7 

4.5 

4.3 

4 .5 

4.1 

2.5 

1 .5 

-.3 

-.2 

1‘.0 

1.0 

- .1 

-.2 


.0 

2.1 

1 .8 

3.2 

3.8 

1 .0 

.2 

- 1.9 

- 2.0 

- 2.1 

- 2..0 

-1 .7 

- 2.0 



.0 

.4 

1.3 


-4.2 

- 2.1 



'.0 

-1 .4 

- .9 


-1 .9 

- 1.9 

- 2.3 

-.1 

-1 .3 










-.9 

- 3.3 

-.9 

3.7 

-1 .9 

1.6 

- 2.3 

- 1.7 

-•9 






.1 

4.5 

?.6 

.3 

8.8 

11.6 

8.5 

2.6 

-.4 






3 .9 

5.4 

5.9 

9.8 

10.5 

5.3 

4n.8 

-.0 

- 3.2 







- 1.5 - 1.9 

-.9 -.6 

-.8 -.3 


a^ 

o 


VELOCltlES lN T DIRECTION 


38 







- 4.3 

.0 

.0 

.0 

37 







- 1.9 

.0 

^0 

.0 

36 







-.5 

.0 

^0 

,.0 

35 







.1 

.0 

>0 

.0 

3 <i • 







.5 

.0 

.0 

^0 

33 







.9 

.0 

.0 

.0 

32 







1 .7 

.0 


'.0 

31 







1 .8 

’.U 


4 U 

30 







2.4 

.0 

;.U 


29 







2.9 

.0 

;.0 

'.a 

28 







3.8 

.8 

1.5 

1.3 

27 






.5 

3.4 

2.2 

2.9 

3.0 

26 






.8 

3.6 

3.7 

4.2 

3.9 

25 






1.7 

3 .8 

5 . 1 

4.1 

3.7 

28 






2.5 

3.7 

5.3 

3.3 

2.5 

23 





.7 

2.3 

3.3 

6.1 

2.9 

5 . 3 - 

22 





1 .9 

2 .8 

3’.0 

4.4 

7.2 

7.3 

21 





2.9 

3.2 

2.6 

8.1 

8.2 

8.5 

20 




1.3 

3.1 

1 .8 

4 .7 

10.9 

9.5 

9.9 

19 




2.9 

1 .8 

4 .1 

4 .8 

12.6 

10.6 

10.1 

18 




2.6 

3.8 

4 .5 

4 .6 

12.8 

11.5 

9.6 

17 




2.5 

4 .7 

5.1 

4.9 

11 .6 

10.6 

9.2 

f £ 

■' 



2 .G 

5..0 

5.7 

5.5 

11.1 

9.8 

8.9 

15 




2.7 

5 .9 

6.5 

6.9 

9.7 

8.6 

8.2 

l«t 




2.9 

5 .7 

6.6 

9.0 

7.3 

7 ,.U 

6.8 

13 




3.2 

5.6 

6.1 

8.4 

5.5 

6.4 

5.9 

12 




3.7 

5 .8 

5 .9 

7.5 

5.1 

5.6 

4.9 

U 



-.8 

2.9 

5.4 

5.5 

7.4 

4.2 

4.7 

4.6 

10 



->8 ' 

3.1 

5.4 

5.4 

7.1 

3.9 

3.5 

4,.0 

9 


J 1 

4 .1 

3.4 

3.3 

3.3 

6.2 

4 .6 

3.3 

2.8 

8 


.0 

-3 .0 

- 3.7 

- 3.3 

- 2.4 

7.8 

4.3 

2.6 

1.9 

7 



‘.0 

- 2.0 

-3 .6 

- 2.4 

7 .7 

2.5 

.8 

.4 

G 




.0 

-4 .4 

-1 .1 

8 .2 

.9 

- 1.4 

- 2.1 

S 





- 4.8 

-1 .2 

2.5 

•0 



4 


.7 

-1 .0 

- 1.9 

-2 . 9 

-1 .2 

3 .0 




3 


1 .4 

- 3.4 

.0: 

.0 

- 9.3 

3.1 

- 5.8 

.3 

- 2.0 

2 


3.7 

1 .7 

-.2 

-3 .6 

.8 

6.7 

- 6. 6 

- 4.6 

- 4.6 

1 


8.2 

.6 

- i.s 

-1 .0 

4 .2 

.8 

- 8.4 

- 10 ,. 0 

- 13.2 


,.0 

2.2 



^0 

4 .2 



>0 

5.8 



..0 

7.0 



*.□ 

7.7 



i .0 

8.1 



*.0 

8.4 



1.0 

8.8 



.0 

9.3 



.0 

9.8 



2.4 

7.1 



2.6 

4.9 



2.9 

2.6 



2.2 

3.7 



4.6 

4.2 



4.7 

4.3 



5.1 

3.1 

!.0 


5.3 

1.4 



5.6 

-.0 



4.9 




6.1 




7.1 

1.6 



7.5 

2.3 



7.1 

3.1 



5.3 

4.2 

.9 
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E. System Verification: 

1. The first step in the verification process is checking the 
tide heights at each tide station with actual data. Actual 
data obtained from tide charts are plotted with the data 
predicted by the model. Comparison of the tidal amplitudes 
and phases between each set of data is then made. An example 
is shown in Figure 29. 

2. The second step in the verification process involves the com- 
parison of discharge rates at each pass within the system. Field 
measurements taken by the Corps of Engineers at Main Pass and 
Cedar Point are plotted with the model predicted values. Figure 
30 and 31 illustrate this method. 

3. The third step in the verification process depends on the suc- 
cessful completion of a cyclic material balance. In this case 
all of the input terms contributing to system mass minus all 
of the mass output terms should be close to zero for the tidal 
cycle (under normal conditions). Input terms include river 
discharge rates and positive discharges at each pass, while 
output terms include negative discharges at each pass and/or 
river or channel. Table summarizes the material balance 
concept for a condition in which bay behavior was normal 
(i.e. no unusual wind, tidal or river flow conditions). 

4. The last step in the verification procedure involves checking 
the flow direction and speed at various locations within the 
system. This method is dependent upon the availability of 
velocity data for point within the bay and, perhaps in some 
cases where surface or bottom data are taken, the interpreta- 
tion of data to conform to the model base (average over depth 
and within the grid area). In most cases, this comparison 
provides trend assessment rather than a specific, quantitative 
check on the system. An illustration is shown in Figure 32 
and Figure 33 for the entire bay. Data used were taken from 
McPhearson (10) and includes ebb flow and flood flow conditions. 

F. The Use of the Hydrodynamic Model: 

✓ 

The hydrodynamic model developed for the Mobile Bay system has been 
used in several studies. Results show that the model has the capability 
of handling a wide variety of problems and of providing important informa- 
tion about the system over several different time frames. A list of 
research studies shown below will give users an insight into several 
model capabilities. 

1. Hydrodynamic paramettic studies in 1972 by Hill and April ( 7 ) 
were conducted to determine the affects of wind direction and 
speed, river flow rate changes, channel effects, convective 
acceleration effects and Coriolis effects on Mobile Bay water 
movement. 



Tide Height (ft.) 


Legend : 


Actual Data 


0““^ — Model Values 



1230 2030 0430 1230 2030 0430 1230 


Military Time 


Figure 29. -- Comparison of Tidal Amplitude and Phases at 
a Location in Mobile Bay. 
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2. Mobile Bay salinity distributions were studied using the 
hydrodynamic model to provide average velocity and dispersion 
data on a tidal cycle average time frame. This study in- 
cluded the investigation of wind effects and river flow rates 
on salinity patterns within the bay.. An extension of this 
work could relate the flushing characteristics of the bay 
during seasonal variations. (7) 

3. Total coliform group distribution studies were studied by 
Liu and April (8 ) in a non conservative species transport 
model. The hydrodynamic model was used to furnish net 
velocities over the tidal cycle and maximum velocities for 
each water cell within the system. Dispersion coefficients 
were also provided in this study. 

4. April and Brett (9 ) used the hydrodynamic model to provide 
current data for isolated sections of the bay in which 
maintenance dredging activities were concentrated. In these 
investigations, transport and distribution of sediment brought 
up from the bay bottom during dredging were predicted using 

a modified conservative species model. This model included 
resuspension of sediment from the bay bottom, settling and 
hindered settling mechanisms. Also studied were trend 
behavior of sediment transport over long periods including 
an assessment of the impact of Mobile ship channel construction 
on bay bottom configuration and sediment deposition patterns. 

5. The hydrodynamic behavior of Little Dauphin Bay has been 
investigated using the hydrodynamic model. The initial goal 
of the project is to predict the distribution and flushing 
characteristics of pollutants within this restricted tidal 
water body. 

In addition to the above on-going investigations, the hydrodynamic model 
will also be utilized in the following future studies; 

1. Storm surge effects on estuary system behavior. 

2. Continued sediment and turbidity distribution studies. 

3. Adaptation of remote sensing data collection methods to 
mathematical modeling techniques. 

4. Mississippi Sound-Mobile Bay interactions at Cedar Point. 

5. Interactive coupling of the hydrodynamic model with the 
various transport models, energy model and other related 
models. 
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G. Limitations of the Hydrodynamic Model: 

At the present time the hydrodynamic model is operational and 
is capable of handling a wide variety of problems. However, there 
are some restrictions which are placed on its formulation to make 
it functional at an economic level. Plans to continue investigations 
to modify the model and thereby remove some of the restrictions are 
underway. These modifications will thus provide a model suitable 
of an even wider range of applications for which bay system behavior 
is needed. 

In its present form the hydrodynamic model is restricted by the 
following; 

1. Only rigid (fixed) boundaries can be handled in the current 
hydrodynamic model. Land areas which could become flooded 
during high water periods are not accounted for (except in 
the river -marsh area in the northern bay) . 

2. The model assumes a constant density throughout the system. 
Thus density induced currents are considered negligible. 
Simulation of salt wedge effects (stratification) is possible 
in those portions of the system where this behavior is known 
to exist, 

3. Wind velocities are restricted to values below 25 knots. 

This condition is partly related to the ability of the model 
to handle only predesignated water cells under restriction 
number 1 above. 
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VII, Appendix A 

This appendix contains a detailed program listing of the 
general hydrodynamic model. 
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♦ •61* f '82' r 'B3' 

f ' 84 ' 

f '05* 

, • 86' 

f ' 87 • f 

•08' ' 

•89' 

• ' 90 ' r 

056 


♦ '91* » '92' » '93* 

1 • 94 ' 

/ ' 95' 

M96' 

f '97' » 

•98' * 

•99' 

, *100*/ 

05o 

c 









056 

c 

PRINTER 









b 


IP = 6 
Tape 

ITaPE2 =2 
I Tapes = 3 
reader 

IMP = 5 
DC = U 

PI =. 4.4ATAM(1. ) 
AMIKAD = laO./PI 
DO 5 I = Ir 20 
lUGRU) = 0 . 
COJlTlNOE 
UaM = 4 
ILP = 3 


MNgHOOIO 

MN 6 H 0 U 20 

MMgHOOSO 

KNGH0040 

MNGH0050 

MNGH0060 

MNGH0070 

MflGMOOSO 

MNGHU090 

MNgHOIOO 

MNgHOIIO 

MNGH0120 

MNGH0130 

MNGH0140 

MN6H0150 

MNGH0160 

MrJGH0170 

MNGH0180 

MNGH0190 

MN6H0200 

MM6H0210 

HNGH0220 

MNGH0230 

MN6H0240 

MNGH025Q 

MNGH0260 

MNGil0270 

MNGH0280 

MNGH029Q 

MN6H0300 

WNGH031D 

MN6H0320 

MNGH0330 

MNGH034Q 

MNGH0350 

KNGH0360 

MN6H037a 

MNGH03fi0 

MNGH0390 

KMGH0400 

WNGH0410 

M(iGH042 0 

MNGII0430 

MNGH0440 

KNGn04bO 

MNGH0460 

MNGHD470 

MNGH0480 

HNGH0490 

MNGH0500 

MNGH051Q 

MNGH0520 

MNGH0530 

MN.G!l0 540 

MrJGH0550 


reproducibility of the 
ORIGINAL PAGE IS POOR 
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OOOOSo 

OOUJb? 

OOUO'jii 

uouubv 

uOuuoU 
oouuul 
OOUCbil 
UOUUbo 
UOOOoH 
OOoUob 
DOUUuu 
l)OUOo7 
UOuUcxJ 
UUUOo-:^ 
00U070 
uOuUTl 
00UU7^ 
OOUU73 
000074 
0 U 0 o 7 o 
U0o07o 
000077 
000070 
uo0u7y 

ooooao 

OOuuBi 

ooooa^ 

ooooao 

000064 

OOOOdO 

UOuCBci 

0Ou087 

ouoOao 

000009 
0OUU90 
0 0u0<31 
000092 
OOOuVO 

U0uU>J4 

U0u09b 

ooouyo 

000097 
00 0090 
UUU099 
OUUiUO 
OUUlOi 
000102 
OOUlOO 
000104 
OOUlOt) 
OUOlOb 
U0O107 
OUOlOo 
oooloy 

ouuiitj 

UOOlll 

OOolld 


05b 

Ubo 

05o 

Ubo 

Ubo 

9 

05rj 

C 

Obb 
Ubo 
05 o 

10 

Obb 

Ubb 

05b 

11 

Obo 

Obo 

U56 

L 

U5b 

Ubo 

30 

ubo 

056 

31 

056 

C 

Ubo 

Obo 

40 

05b 

05b 

41 

05o 

05 d 

U5fa 


056 

056 

05o 

056 

C 

05fa 

05b 

056 

6u 

05b 


05b 

U5o 

61 

056 

C 

056 

05o 

Ob 

05b 

C 

05b 

056 

056 

05b 

C 

Obo 

Obb 

74 

U5b 

C 

U5b 

Obb 

05b 

80 

Uoo 

Obb 

Obo 

B1 

056 

C 

Ubb 

05o 

90 


1SE.T1 r 1 
lbt,l'2 = 1 
lbt.T3 = 1 


V.KI rE(lP,9) 

FOHriATi * 1 • lOf * ii»puT Data check * 
keau It. ppooram cohtrol card 

REMD(IttP.lO) IFMT, ICORF, ICONVr lORIGr 
ilrtODFe ICKr IMMDPf INC 
FORMAT! 1112) 

Wi<UE(lP»U) IEMT» ICOHF. ICONV, lORiG. 
IIMODF. ICKf IKriDFr INC 
FORMArc •.1112) 

IF (IFf'T.EO.O) IFMT = 6 
IFh,.U 7) r I5 PaCE{IFMT-4) 

REaU id geometrical information 
RE rtuai.p.ao) iLB. ipow. IV, lu. ir, nc, 
FORMATi7I5) 

wRlTEtIP.31) IL-P- IROWi IV, IU, IR, NC# 
format (• ..TIS) 
read id physical information 
R£AD(H-P t40) W, THETa. DELS, DELT , TIM, 
FORMAT (7F10. 3) 

AlUTEt IP,«i) W, THETA, DELS, DELT, TIM, 


*** ,///) 

IRMil, ISPUK, IDSCH, 
IRMB, ISPUK, IDSCH, 

NCPR 

NCPR 

SLaDJ 

SLaDJ 


format (t <,'7F10.3) 
FAA = DELT/4. 

FAU = celt /DELS 
FAC = lo. l«l)ELT/DELS 


TIM = TIM»3f>00. 
INC = INC43fc00 


IN A ROW 


READ ID left AHD RIGHT BOUNDARY CELLS OF WATER CELL SEGMENTS 
REAO( INR'bO) ( IBNDL(I) , 1=1, ISB) 

READ! It:P,60) (IbNDR(I) f I = 1.»ISB) 

FORMAT! (16IS)) 

WKlIE!lP,faI) iluMOLd) , 1 = 1, ISB) 

.VRl rElIP,6I) nuNDRlI) , 1 = 1, ISB) 

FORMAT!' ',1615) 

lRc.P!l) H OF PRINT REPETITIONS FOR EACH SEGMENT 
DO 65 I = 1,1 SB 

lRtP!l) = inNOB!!) - IBNDLII) +1 

NUM!I) IS THE SPACE INDENTED BEFORE PRINTING BEGINS IN EACH RO.t 
NUi'1!l) = lIPNDLCl) -1 )*TFMT + 5 

00 74 J =2, ISP 

HR = IBNDLIJ) - IUl'fDLlJ)/NCPR*DCPR - 1 
PRINT. rO,!MAT ISF!IFMT).I, THEREFORE KK^IFMT 
NiJM( J) = Kt<» IFt^T + 5 

REaD^ID^CELLS THAT REQUIRE 7ER0TMG OF CALCULATED VEL. 

KEaHI iDPr IU)) ! II.'C 1 1 j , 1 = 1, IV) 

RLADUDP.ao) IICDII), 1 = 1, lU) 


FORMAT !1 dI 5) 

WRITE! IP.Bi) t IF'CI I J I 1 = 1, IV) 

V.'R'iTE IIP, till Ui:D(I), 1 = 1, K'l 
FOR.'iAl ! ' • ,1615) 

READ ID UlSCHAl'GE RATES Of RIVERS 
kEADHDP,'?0)!Rf(I)> i = l,IRl 
FOilMAl! !!iFlG.b) ) 


MNGI 10560 

HNGH0b70 

MtJGHOSaO 

MUGHObSO 

KNGH0600 

|.)N&llC610 

MNGHD620 

MIIGH0630 

MNGI 1064 0 

PNOH0650 

MNGH0660 

M1JOH0670 

MNGH0680 

MNGH0690 

MNGH0700 

Mr)GH07i0 

MN6H0720 

MNGHBT30 

MNGH0740 

MN6H0750 

MNG! 10760 

MNGH0770 

KN6H07BO 

KNg )!0790 

MNgHOSOO 

MNGH08T0 

MI\IGH0820 

MN6H0830 

KNGHOaLO 

KNGH0850 

MNGH0S60 

MNGH0870 

mNghobso 

MNOH0890 

HNGU0900 

MN6H0910 

MNGt 10920 

MNgHU930 

MNGH0940 

MNGHOgbO 

KNGH0960 

MNGH0970 

MNGH0930 

HNGH0990 

MNgHIOOO 

MNGHIOIO 

MNGH1020 

MMSH1030 

MN61I104 0 

MNGM1050 

MNGH1060 

MNGH1070 

MNGHIOOO 

KlJ&tll090 

HNGHUOO 

MNGHillO 

MN6H1120 


82 


OOU113 

056 


V.RirL(IP.91) CRb(I) . I = 1*IH) 

U0U114 

obo 

91 

FORMaTC <,(0fl0.5)) 

UOullb 

Obb 


00 r(5 I = IfNC 

OOOllo 

Ubb 


AMuY(I) =0.0 

OUOliV 

0 5b 

8b 

2U) = 0.0 

OOullb 

Obo 


DO 9b K = lilSb 

UUU119 

Obb 


KF = IBNDL(K) 

oooiao 

U5b 


KL = U'NDRCK) 

000121 

Obb 

C 

RlaD in tVOiNUIG COtFFICIENTS 

000122 

056 


REAU(IMP» 100) (aMNY(I) .I=KF fKL) 

000123 

05b 

100 

FORMAT! (26F3.Q) > 

00ul24 

05b 

9b 

V,R I TE ( 1 P • 10 1 ) ( aMn Y C I ) * I=KF f KL > 

00ul2S 

Obb 

101 

FORMAT!* *,!26F3.0)) 

00012b 

05o 


Do llu I = l.NC 

0UU127 

Obb 

no 

AHNYCl) = AMMY!1)/1000. 

000128 

Obb 


DO 120 K = 1»1SB 

000129 

Obb 


Kl- = IBNI)L(K) 

000130 

Obo 


KL = IuNDR!K) 

000131 

Obb 

C 

RLAD IN DEPTHS OF ,/ATEP CELLS 

000132 

05b 


READI1!;P*140) 121 I) . IrKF.KL) 

OOU133 

05o 

140 

format 1 12OF3-.0) ) 

DUU134 

U5o 

120 

WRITEUP, 141) !4(I) . 1=KF,KL) 

0OU13b 

056 

141 

FORiIaT!' *i!26F3.0)) 

U0ul3o 

05o 


DO 15b K = 1, Isa 

000137 

05d 


IbTRT = lONDLlK) 

000138 

056 


IJUIT = IHmDRIK) 

U0O130 

Obb 


DO 154 I - ISTRT, IQUIT 

ooomo 

056 

154 

211) = Zin+SLADJ 

00 0141 

05b 

155 

CONTINUE 

000142 

05b 

C 


000143 

0 5b 

C 

initiaeization 

000144 

Obo 

C 


00014b 

05d 


DO 170 I = 1*NC 

00ul4o 

U5o , 


UU) = 0.0 

000147 

Obo 


V!l) = 0.0 

000148 

056 


UPtI) =0.0 

00U149 

05b 


VPll) = 0.0 

OOUiSU 

05o 


till) = 0.0 

000151 

05b 


D(I) = 2!D 

00U152 

056 


VMAGII) =0.0 

000153 

056 


BC.TAU) = 0.0 

000154 

056 


SUMX ! I) = 0.0 

00015b 

Obb 


SUMY m = 0.0 

OUOlSu 

Obo 


SUMXNTII) =0.0 

000157 

Ubb 


SUMYNTd) = 0.0 

000153 

05b 


FX! I) = 1.0 

000159 

05b 


FY(I) = 1.0 

00 0 160 

U5b 

170 

CONtliJUE 


UOUiol U'jo R = 0.0 

000162 Obb C VELOCITIES ''‘OOlFICATIOtl factor 

OOOlbJ OSo IP (1MODF.CO.O) GO TO 190 

000104 Ubb C INPUT NO OF CELLS TO RE MOQl FIFO 

OUUlfjb .O'jo 11 l,\D( 1NR*172) NrfOoFXt NROOFY 

OOulbb U5b 172 F0KMAT(215) 

UOoifo? Obu WrAlTr{lp,i73) NROOFX, NMODFY 

OOUlbb Ubu 173 FORMAT (•' »2I5) 

000169 05u C INHUT CELL rio Velocity MODIFICATION FACTOR 


RNGH1130 

MNGH1140 

MMGHlibO 

MNGH1160 

MNGH1170 

mNgtuioo 

MUGH1190 
MNGH1200 
MNGH1210 
MNGH1220 
HN6M1230 
MMGH1240 
MtJGH12S0 
KNGH1260 
MN6H1270 
MNGH1280 
MNGH1290 
MNGHlsao 
MNGII1310 
MN6H1320 
MN6H1330 
MNGH1340 
NN61U3S0 
MNGl 11360 
KMGH137D 
MNOHlCi'O 
MNGH1390 
MNGHmOQ 
MMGHimO 
MN6IU420 
MNGH1430 
MNG1T144 0 
MNGH1450 
MNGHlRbO 
MNGH1470 
MNGI 11480 
MMGH1490 
HNGHlbOO 

mnguisio 

KNGH1520 

MNGI11530 

MNGH1S40 

MNGIIlbbO 

MNQHlbbO 

MNGH157Q 

MNGH1580 

MHGHlb90 

MNGUibOn 

MNGII1610 

MW6H1620 

HNGH1630 

MNGH1640 

MNGHlbbD 

MNGII1660 

mNgh1670 

MNGH16BQ 

MNGH1690 



83 


000170 

0 5o 

C 

X-PIKECTION 

U0U171 

058 


IF (t.VlODFX.tO.O) GO TO 100 

000172 

056 


HwAn(liJr>.175) (KC(I),VHl).I = ltl4M0DFX) 

000173 

056 

175 

F6KI-lAT(6Ub.F5.?) ) 

OU0174 

U5o 


.-•■rvl ;r;7.P,176) lKC(I),VFtI»,.I = l,fJMODFX) 

000176 

05o 

176 

F3-i • .«( I5fF5.2) ) 

00ol7o 

Ubu 


uo ; 7, I = l.MVOUFX 

00oir7 

05o 


J = KC< I ) 

000178 

05o 

179 

Fa(J) = VFCI) 

000179 

056 

C 

Y-oIHOCTIOii 

000180 

U5u 

loO 

IF (WMODFY.EO.O) GO TO 190 

000161 

05q 


l<8ALl(IiJP,16r3) (KC(I) .VF(I)*Irl,NMODFY) 

000182 

Obo 

105 

F3HRAT (B(I5»F5.2) ) 

00U16.'> 

05o 


W:UTE(IP.10b) (KC<I) ,VF(I)fI=l,NMODFY) 

000184 

05u 

166 

format ( • I r8{ I6»F5.2) > 

OUol8b 

Obo 


00 189 1 = l»NN.ODFY 

oouiao 

05o 


J = KC(I) 

000137 

05o 

189 

FY(J) = VFII) 

OUOlOO 

056 

190 

CoNTlMUE 

OOU189 

050 

C 


000190 

05d 

C 

CALCULATE WIND STRESS 

000191 

0 5o 

C 


000192 

Obo 


IF (W.LT.23.64) go TO 220 

000193 

056 


AX = (.00011 .00025*(i.-23.6446/wn/100. 

000194 

05o 


GO TO 230 

000195 

056 

220 

AK = .0000011 

00ul96 

05o 

230 

CONTINUE 

00U197 

05o 

e 


000198 

Obo 

C 

REKEMbER TO MEt-'TIOiJ Tt(E DIRECTION OF CALCULATION 

000199 

05o 

c 

theta IS MEASUKEO FROM X-aXIS COUNTER-CLOCKWISE 

00O20U 

056 

c 


000201 

05o 


X = AK»W*W»C0S(THETA) 

000202 

056 


Y = AK*W»W*SIN(THETA) 

0 00203 

05o 

c 


000204 

056 

c 


000205 

058 

c 

Calculate CORIOLIS force 

000206 

056 

c 


0002J7 

Obo 


COR = 0.0000727 

00O2U8 

OSo 


IF (ICORF.EO.O) C0R=0.0 

UUu209 

Obo 

c 


00O210 

Obo 

c 

PRINTING CONTROL FOR MANNING FRICTION FACTOR AND INITIAL DEPTH 

U0O211 

05o 

c 


000212 

056 

c 

■ ' * 

00u2l3 

056 


IF ( IMNDP.EQ.O) - GO TO 245 

000214 

0 So 


call PRINT1(4»AMNY» 

000215 

Obo 


IM'IANNIUG FRICTION FACTOR 

000218 

056 


CALL PRlNTl(?..2r 

000217 

Obo 


I’JLPTH OF water FrON DOTTOM OF BAY TO REFERENCE PLANE 

000218 

Cbo 

245 

06 2SU I = l.MC,, . 

000219 

Obo 

250 

AMNY(I) = 14.57*(A,rtNYCI) *»2) 

000220 

Obo 


nC : INC 4 T 

000221 

05o 


Call 8DUND1 i INPi IP. TMlLi. ITSN.Tp ,HDB» FAB .FAC . DELT . DELS. NCPR'. R ) 

000822 

Obo 

C 

READ TAPE 

000223 

u5o 


IF . i ICK.EO.O) 00 TO 450 

000224 

OSo 


P.n.lNJ ITAPC3 

000225 

05o 


READ(ITAPE3.400)(H11). I=1.NC) 

000228 

05i. 


RLADHTAPE3.400) (U(I)f 1=1. NC) 


MUGH1700 

MNGH1710 

MNGH1720 

NJGH1730 

MNGH1740 

MI'IgH1750 

MNGH1760 

MN6H1770 

MNGH1780 

MNGH1790 

MNGH1800 

MMGIU810 

KNGH1820 

MNGH1630 

MI'IgH1640 

MNGH18S0 

MNGM1660 

MNGH1870 

MN&H188C 

l/iNGHiagO 

SNGHlgOO 

f-iNGHigiD 
MNGH1920 
HN6H1930 
MNGH1940 
MNGH1950 
KNGHI960 
MNgH1970 
RNGHigSO 
RN6H1990 
MN6H2000 
KNGH2010 
MNGH2020 
ntlGH20 30 
MNGH2040 
MNGH2050 
MNGH2060 
MN6H2070 
MNGH2080 
MNGH2090 
KNGH21Q0 
MNGH2110 
MNGH2120 
MNGH2130 
MNGH2140 
' )MNGH2150 
MNGM2160 
* )MNGH2170 
MNGH2iaO 
RN6H2190 
RNGH2200 
HNGH2210 
MNGH2220 
t-'NGt 12230 
MNGH2240 
MMGH2250 
- MMGH22.60 


REPRODUClBrLrrY OP TBM 

ORIGINAL PAG!?] IS POOS 
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00U227 

00U229 
U0O230 
000231 
000232 
000233 
U0023H 
00o23b 
000230 
00C237 
00023d 
000239 
000240 
000241 
000242 
000243 
000244 
00024b 
00O24& 
000247 
00o24d 
000249 
00O250 
000251 
OOU252 
000253 
000254 
000255 
000256 
000257 
000256 
000259 
o0o260 
000261 
000262 
0002o3 
000264 
00U265 
000266 
00 0267 
000266 
00U269 
00 0270 
000271 
000272 
000273 
000274 
000275 
000276 
000277 
OOU276 
000279 
OOU260 
000261 
000262 
000263 


056 


(Jbo 

400 

U5b 

450 

05b 


056 


Ubb 


05o 


056 


056 

460 

056 . 


056 

470 

05o 


056 


056 


056 


056 

460 

OSo 

SCO 

056 


056 


056 


05b 


05b 


OSo 

54 0 

056 

550 

056 


066 


U5o 


056 


05b 


05b 


056 


05b 


056 


05b 


OSo 

590 

056 

C 

05b 


05b 


05b 


05b 


05b 


05b 


056 


05b 

oUO 

05b 

c 

056 


05b 


U5o 


056 


056 


05b 

610 

0‘jU 

615 

Obo 


OSb 


05b 


05b 


OSo 



Hh/O(lTAC>E3.40tJ)(V(l)» 1 = 1. MC) 

FOH(.',aT(T2,1HF7.3/T2. 11F7.3) 

COlJTlrJbE 

IF (I5E;n.tG,2) GO TO 460 . « 

CAl.L \/t Ll (INP. IF'. Ii<f‘'Bfl^CPH.riC.Df:LT.D!:LS,I0Kl6.1C0NV.AMtJY.C0R.IU. 
llv. I5D. I6C.15tl,FX,FY,FAA.FA6.FAC,X. Y) 

I5LT1 =2 
GO TO 470 

CALL VF.L2(IIJP. IP« liWB.IJCPN.f'IC.OELT, DELS, IORIG.ICONV.AMNY. COR, lU. 
lIV,IUU,I(iC.I56.FX,FY.FAA,FAB.FAC.X,Y) 

IF (ISPDK.cO.O) GO TO 500 

IF (I5ET2.EG.2) 60 TO 480 

C/.LL btiAPKI ( Itif’, IP ^NCPR .DELT .DELS ,FAA .FaC .COR , X . Y) 

ISLT2 =2 
GO TO 500 

CALL ShANK2(INP, IP.IJCPR.DLLT.DELS.FAA.FAC.COR.X.Y) 

continue 

IF (IDSCH.EO.O) GO TO 550 

IF (ISET3.E0.2) GO TO 540 

CALL UISCHKINP, IP .DELS . DELT. VOL, DlSR .SDISRP. 5DISRN. LDSCH) 

IbtT3 =2 

GO TO 550 ^ . ■ . , 

CALL UISCH2(INP. IP, DELS, DELT, VOL, DISR .SOISRP, SUISRN, LDSCH) 

COIiTIMUE 

DO 590 J = 1 rise 

ISTRT = IBiJOLiJ) 

lUUIT = I0 iNDR(J) 

DO 590 I = ISTRT, IQUIT 


XaX = SUl^Xd) 
xxp = ABSfUPtI) ) 

IF tXXP.GT.XXX) SUMX(I) = XXP 
YTY = SUHY (I) 

YYP = ADS( VP(I) I 

IF (YyP.GT.YYY) SUMYTIl = YYP 

CONTIiJUE 

CONVERT VP TO V AND UP TO U 
DO 600 J = l.ISB 
ISTRT = iBNnL(J) 

UUIT = IBUDRIJ) 

DO bUO I = ISTRT, lUUlT 
SUWXNT(I) = SUN'.XNTd) 4 UPd>/D(I> 
SUMYHTd) = SUMYNTd) + VPd)/D(I) 
udi = uPd) 

Vd) = VPtI) 

CALCULATE TIDE HEIGHT IN BAY 


DO 610 J = l.IbB 

IF ( luUDL( J) .GE. (NC-NCPR) ) GO TO 615 
IbTRT = lUNDL(.J) 
lOUIT = IBNDR(O) 

DO 610 I = ISTRT, IQUIT 

Hd) = Hd> 4 FAB*{UP(n+VP(r)-UPd+l>-VPd4NCPR) )+R»DELT 


continue 

IF* (Ai-lODdMIL,1440. ) .LT, 0.0001) TMiLrO.O 

T r T + OEUT 

TMIL = TF1IL + LELT/60. 

Call bOUFJD2(INP,IP, TMIL.ITSN.Td.HdB, fab fFAC.UELT. dels, HCPR.R) 

n = T 


WNGH2270 
MNGH22H0 
MNGH2290 
MNGH2300 
MNGH2310 
MNGH2320 
MN6H2330 
MNGH2340 
MNGH2350 
F’NGI 12360 
MNGH2370 
MNGH2380 
MNGH2390 
MNGH2400 
MN6H2410 
MNGH2420 
MN.GH2430 
(.•.NGH2440 
KNGH2450 
MNGH2460 
MN6H2470 
MNGH2480 
MNGH2490 
MN6H2500 
MNgH 2510 
MNGH2520 
MNGH2530 
MNGH2540 
MNGH2550 
MNG! 12560 
MNGK2570 
MNGH2580 
MNGH2590 
MNGH2600 
MNGH2610 
MNGH2620 
MNGH2630 
KNGH264a 
MN6H2650 
MNGH2660 
MNGH2670 
MN6H2600 
F-.NGH2690 
mNgH2700 
MNGH2710 
F1NgH2720 
MNGH2730 
mNgH2740 
MNGH2750 
MNGH2760 
MNGI 12770 
F,MgH 2780 
.KNGH2790 
MN6H2800 
W'lGH281 0 
FNGH202O 
MNGH2830 
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IJELT.L SA.’fl.'jUUl 

EU0T7 «LU170 *5: ba- (24 . ) 


OOUUUi 
U0UU02 
U0(;00 J 
OOOOU4 
OUOUUb 
OOOO'Uo 
0OU007 
OOOOUo 
000009 
UOUUIO 
OOOOII 
U0UU12 
00UU13 
000014 
000015 
OOOQlo 
U0U017 
UOUOIB 
000019 
000020 
000021 
000022 
00UU23 
000024 
000025 
000026 
000027 
OOU02B 
000029 
000030 
OOU031 
000052 
000033 
000034 
000035 
000030 
000037 


024 

024 

024 

024 

024 

020 

020 

020 

02U 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 


b 

10 


51 


52 

54 


VEHSI0IM.2B.7b 

♦ 

50UR0UTINE PRIliTldXYiSTOREl.X) 

DIMENSION 5TORt*lC2500) f 5TORE2(2500) * X(ll) 

COMMON /A/ Z(2b00)» H(2500). D(2500)i lBNDL(200>r lBNDR(200)i 
INUM(IOO) i IRE(M200> 

COMMON /C/ NCPI(» NCr IROWt ISB» INP» IPr 1PRNT(6)* lFRM(10)t 
UFLD(4), INDENT (100) 

IPASS =1 
60 TO 5 

ENTRY PRINT2(IXYfST0REl»ST0RE2iX) 

IPASS =2 
WRITE(IP.IO) X 
F0HMAT(il*.T2rllA6) 

IFRM(3) = IPRNTC5) 

IFR(J|(0) = IFLDIIXY) 

JK = ISh 
M = 0 
DO 70 K = 

M : IROW ' 

IFKM(2) = IPRNT(l) 

JL = JK 
JL = JL-1 

IF (JL.EQ.O) 60 TO 52 
IF ( ( lBfJDL(JL)/NCPR+l) .EQ.M) 60 TO 01 
JF = JL 
JL = JL + 1 
Nil = NUM(JL) 

IFHM(4) - INDENT(NM) 
lU = IREP(JL) 

IFUM16) = INDENT(IR) 

ISTRT = IBNOL(JL) 
lOUIT = 1BNDR(0U 
IF (IPASS. EQ. 3) GO TO 55 


l,IROW 
• K + 1 


000038 

020 


WRITE! IP. IFRM) HdSTOHEll 

000039 

020 


IF (JL.EG.I) GO TO 56 

000040 

020 


IF (JL.EQ.OK) GO TO 58 

000041 

020 


GO TO 56 

00UU42 

020 

55 

wRITEdPf IFRM) Mf(5T0RE2 

000043 

020 


IF ! JL.EQ. 1) 60 TO 70 - 

000044 

020 


IF (JL.EQ.JK) GO TO 59 

000045 

020 

56 

IFUM(2) = IPRNT(4) 

000046 

020 


60 TO 54 

000047 

020 

58 

IF (IPASS, EO.l) GO TO 60 

00004B 

020 


1FRM(2) = IPRNTd) 

000049 

020 


IPASS =3 

000050 

020 


JL = JF 

000051 

020 


GO TO 54 

0UOO52 

020 

59 

IPASS = 2 

0OU053 

020 

60 

JK = JF 

OOU054 

020 


IFRM(2) = IPRNTd) 

000055 

020 

70 

CONTINUE 

00UO56 

020 


return 

000057 

02U 


END 


IPRTOOIO 

1PRT0020 

1PRT0030 

1PRT0040 

1PRT0050 

IPRT0060 

1PRT0070 

IPRTOOBO 

1PRTD090 

IPRTOIOO 

IPRTOllO 

1PRT0120 

1PRT0130 

1PRT0140 

1PRT0150 

1PRT0160 

1PRT0170 

IPRTOISO 

1PRT0190 

1PRT0200 

1PRT0210 

1PRT0220 

1PRT0230 

1PRT0240 

1PRT0250 

1PRT0260 

1PRT027O 

1PRT02BO 

1PRT0290 

1PRT030D 

IPRT0310 

1PRT0320 

1PRT0330 

1PRT0340 

1PRT0350 

1PRT0360 

1PRT0370 

1PRTD380 

1PRT0390 

1PRT0400 

1PRT0410 

1PRT0420 

1PRT0430 

1PRT0440 

1PRT0450 

1PRT0460 

1PRT0470 

1PRT04S0 

1PRT0490 

1PRT0500 

lPRT05)0 

1PRT0520 

1PRT0530 

1PRT0540 

1PRT0550 

1PR(U560 

lP«T0570 


EI4U LLT. 
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«E.!-T»L ^Arl.S>UEl2 

tLlUT? i<Lib7U U1/29-17;‘»6:02-<52. ) 

OOOUOl 052 C 

00UU02 052 C ♦ ‘ . 

OOUOOS 052 C VfRSlOIJ 1,28.7b 

OOU0U4 052 C ♦ 

000005 052 C 

OOOOOo 046 SUbROUTIME REPRTl ( IP, IIJP, ILP, LCSCH, inSCH» IRf Kli. I ^AM» IT, I TC» 

OOU007 04o 1I!JC,DELT,0ELS,T,T1 ,vi,V 0L,HC, ITAPE2, lTAPE3,$,i,1.) 

000006 048 OlhENSION VOL(20) . Rfi(20), RF(20) 

000009 048 common /A/ Z(2800), Ht2500), D(25oo>, IBNDU200), IBNDK<200) 

OOuOlO 046 INUMUOO), IREP(2001 

OCUUH 04a COMMOfJ /»/ VMAO(2500), RETA(2500), AHK2500), AH2(2500), 

000012 046 1SU(.',X(2500), 5UMY(2SOO), Sul-IXNT ( 2500 1 , SUMYNT ( 2500 ), U ( 2500 ) , 

00u0l3 048 2V(^500) ,UP(2500) ,VP(2500) 

000014 048 IP (ILP ’T.3) 00 TO 90 

00U015 048 IP tIOLUi.EO.O) GO TO 30 

OOOOlo 046 DO 10 I = l.LDSCH 

000017 040 10 WimE(lP,2U) I, VOU I ) 

OOoOia 048 20 format (T 2 , 'VOLUME FLOW RATE AT LOCATION (•,12,*) =*,F15,2, 

000019 048 12.x, *00. FT. /SEC ) 

UOU020 048 JO continue 

000021 043 UAM = UaM + 1 

000022 048 IF (IT.EQ.ITC) GO TO 35 

000023 048 RETURN20 

000024 048 35 CONTINUE 

000025 048 DO 40 1 = 1,IR 

000026 048 Rt- = R3(I)*DEL5 

000027 046 WRITEdP.Sd) 

000028 048 3o FORMAT(* •,/) 

000029 048 40 WHirE(IP,50) I, RF'lI) 

000030 048 So FORMAT ( T2 , • R I VER - • , 12, 2X , • FLOW RATE =',F10.2,2X, 

000031 048 1*CU.FT./5EC* ) 

0OU032 048 ITC = ITC + INC 

0UU033 046 CALL PRINT1(3<H, 

000034 046 I'ilEIGHT FROM REFERENCE PLANE TO SURFACE 

0OU035 048 Cai_L PRINT1(2,U, 

OOUOJo 048 1* VELOCITIES IN X DIRECTION 

000037 046 CALL PRiNTl(2rV, 

OOU038 048 1' VELOCITIES IN Y DIRECTION 

000039 04o CmLL PRINT2 ( 2 , VMAG, BETA, 

000040 048 2 'MAGNITUDE AND DIRECTION OF VELOCITIES 

000041 048 CALL PR INT2 ( 4 , AHl , AH2, 

000042 048 I'CONVECnVE ACCELERATION 

OOU043 048 90 CONTINUE 

000044 046 ENTRY REPRT2 ( IP , INP. ILP .LUSCH, IDSCHdR • RfT, I JAM, 

000045 048 11T,1TC,INC,DELT,DELS,T.TIM,V0L,NC,1TAPE2,1TAPE3,S,4,3,) 

000046 048 DO 95 I = 1,NC 

O0U047 048 9L D(I) = Z(l) + h(I) 

000048 048 IF (ILP.LT.3) RETURN21 


000049 

048 


IF (T.LT.TIM) RETURN22 ’ 

000050 

048 


DO 100 I := 1,NC 

000051 

048 


SUMXNT(T) = SUPXNT(I)»nELT/TIM 

000052 

048 

100 

SUMYNT(I) = SUi;yNTa)*OELT/TIM 

000053 

046 


CALL PHI(JT2(1,SUMX,SUMY, 

000U54 

048 


I'UISPEHSION COEFFICIENTS 

000.Q55 

048 


Rc-WIND ITaPE2 

000056 

048 


RE.vINO ITAPE3 

00 0057 

046 


WtUTE(ITAPE2,110) (SUMX(Il,I=l,flC) 

000058 

048 


WRITE [ ITAPE2 * 110 ) ( SUMY (I ) , Isl, NC ) 

000059 

048 


v.R1 TE( ITAPr.2,110) (SUMXNT(I) ,I=1,NC) 

OouudO 

048 


WKITE.( ITaPE2,110) (SUmYNTC I ) , I = 1,NC) 

oouroi 

048 

UO 

FOKMATl (T2f7F8.0,/) ) 

000062 

048 


WRITE(rTAPE3,115) (H( I) ,1=1 ,NC) 

00uU63 

048 


WRITEdTAPES, 115) (U ( I ) . 1=1 , NG ) 

000064 

046 


WKTTEITTAPES.IIS) (V(I),I=1,NC) 

OOuOOb 

048 

1 15 

F0RMATIT2, 10F7.3/T2,11F7.3) 

OOO066 

048 


DO 117 I = 1,NC 

000067 

048 


bUMXNT(I) = 5U:.,XHTtI)#1000 

OOU060 

048 

117 

SUMYNTtl) = SUt.YUT(l)*1000 

00 g 069 

046 


CALL PRlNT2a,SUMXNT,SUMYNT, 

000070 

048 


2<NET velocities IN X-Y DIRECTION ♦1000 

ei'0071 

046 


DO 120 I = IrNC 

UU0072 

048 


SUMXNTtn = SU.MXNT(I)/1000. 

000073 

048 

120 

SUMYNT(T) = SUPVtJT(I)/1000. 

00U074 

048 


RETURN21 : 

000075 

048 


END 


2RPT0010 
2RPT0020 
2RPT0030 
2RPT0040 
2RPT0050 
2RPT0060 
2RPT0070 
2RPT0080 
2RPT0090 
2RPTD100 
2RPT0110 
2RPT0120 
2RPT0130 
2RPT0140 
2RPT0150 
2RPT0160 
2RPT0170 
2RPT01B0 
2RPT019O 
2RPT0200 
2RPT0210 
2RPT022D 
2RPT023Q 
2RPT0240 
2RPT0250 
2RPT0260 
2RPTU27D 
2RPT0280 
2RPT0290 
2RPT0300 
2RPT0310 
2RPT0320 
2RPT0330 
' )2RPT0340 
2RPT0350 
• )2RPT0360 
2RPT0370 
« )2RPT0360 
2RPT0390 
' )2RPT0400 
2RPTC410 
• )2RPT0420 
2RPTD430 
2RPT0440 
2RPT0450 
2RPT0460 
2RPTG470 
2RPT0400 
2RPTP490 
2RPT0500 
2RPT0510 
2RPT0520 
2RPT0530 
' )2RPT0540 
2RPTD55Q 
2RPT0560 
2RPT0570 
2RPTUS80 
2RPTf)590 
2RPT0600 
2RPTU610 
2RPT0620 
2RPTG630 
2RPT0640 
2RPT0650 
2RPTU660 
2RPT0670 
2RPT06BQ 
2RPT0690 
• )2RPT070Q 
2RPT0710 
2RPT0720 
2RPT0730 
2RPT0740 
2RPT075D 


END ELY. 


88 


CJOUUOl 

006 

c 


0UUU02 

006 

c 

'4 

OOUO03 

UUu 

c 

VLKSIOIj 1.28.76 

000004 

00b 

c 

♦ 

OOOOOb 

UOu 

c 


OOUOOb 

005 


subroutine VELKIIlP.IP, IRI-i'J . IICPR i t4C .OELT« DELS. lORIG. ICONV. AMNY 

00UUU7 

005 


ICOK.IU.IV.lnD.lHC.lSD.FX.Kr.FAA.FAl'.FAC.X.Y) 

ouuouo 

005 


dimension AMNY(2500), IHVi 38(2)* rKF!-T(2). FX(2500), 

000009 

005 


1FY(25U0), U!C(250), IbU(25U) 

oooolo 

005 


COMMON /A/ 2(2500), H(2500), D ( 2500 ). lOMOL 1200 ) . I3NDR(200), 

ououu 

U05 


INUt^idOt.') . IFEP(200) 

00U012 

005 


COMMON /D/ VMAG(2500). BETA(2500). AMK2500). AH2(2500), 

000013 

005 


1SUMX(2500) . SUMY (2500) , 5UMXNT( 2500 ) , SUMYNT (2500) , U(2500) » 

00O014 

005 


2V(2500), UP(25(IO)r VP(2500) 

OOuOlO 

005 

c 


000016 

005 

c 


000017 

005 

c 

CALCULATE VELOCITIES IN X-Y DIPECTIOMS 

oouoia 

005 

c 

THIS section is DESIGNED FOR TIC RIVER MARSH BOUNDARY IN (■lOnlL 

000019 

005 

c 

BAY SYSTEM OR SIMILAR 'SYSTEM 

000020 

005 


IFROW =99999 

00OU21 

005 


IF (It(f'B.Eu.O) GO TO 100 

000022 

005 


REA0(IKP.50) (1RMBi3(I).1=1,2) 

000023 

005 

50 

F0R[>1aT(2I5) 

UOU024 

005 


hRITEdP.Si) (IRMt3Bd).I=l,2) 

000025 

005 

51 

F0RMAT( • • .215) 

OUU02U 

005 


RtADdNP,50) (IR(-luTd) ,1=1.2) 

000027 

U05 


READdMP.60) ZRB 

000028 

005 

60 

FORMaT(F10.5) 

000029 

005 


li^'KlTEdP.Sl) (lRMBTd),I = 1.21 

000030 

005 


Wi<ITE(IP,61) ZRB 

00u03i 

005 

61 

FORMAr(» *,F10.5) 

000032 

005 

80 

DU 90 I = 1>NC 

000033 

005 

90 

Dd) =H(I) 4 2d) 

000034 

005 


IFROW = 1RM[TB(1)/NCPR + 1 

00O035 

0 05 


entr y VEL2 ( I NP . IP , IRMB r NCPR , |JC , BELT .DELS . I OKIG . ICiONV . AMNY , COR , 

000036 

205 


llV, IBD, IDC. ISB,FX,FY,FAA>FAC-.FAC4X, Y) 

00003? 

005 

100 

DU 310 J = l.lSf) 

000038 

005 


IF dt)NDL(J).GE.(NC-NCPR)) GO TO 350 

000039 

005 


ISTRT = IBNDL(J) 

000040 

005 


IQUIT = IBNDR(U) 

000041 

005 


UCONVa =0.0 

000042 

005 


VCOMVA =0.0 

00OU43 

005 


DO 310 I = ISTRT. IQUIT 

00UU44 

005 


Ii( = I + NCPR 

000045 

005 


51 = Vd) + Vd + 1) + V(IH) + V(IH+1) 

00004O 

005 


S2 = U(I) 4 U(l4l) + UdH) 4 U(Itl4l) 

00U047 

005 


IF ( dOUIt-ISTRT) .EQ.O) 52 = 0, 

000048 

005 

CCC 


000049 

005 

C 


000050 

005 

c 


OUOObl 

005 


IF (ICONV.EO.O) GO TO 150 

000052 

005 


IF ( (I5TRT/NCPIO.LT.1) GO TO 150 

000053 

005 


IF ( J.0E.1FP.0W) GO TO 150 

000054 

005 


FACC = nELT/(DLLS»(Zd)+Z(l4l> ) > 

000055 

005 


FACD = DELT/(DLLS»(2(I)+2(IH))) 


3VEL0010 

3VEL0020 

3VEL0030 

3VELO040 

3VEL0050 

3VEL0060 

3VEL0070 

3VEL0080 

3VEL0090 

3VEL01UO 

3VEL0110 

3VEL0120 

3VEL0130 

3VEL0li40 

3VEL0150 

3VEL016Q 

3VEL0i70 

3VEL0180 

3VEL0190 

3VEL0200 

3VEL0210 

3VEL022D 

3VEL0230 

3VEL0240 

3VEL0250 

3VEL0260 

3VEL0270 

3VEL0280 

3VEL0290 

3VCU0300 

3VEL0310 

3VEL0320 

3VEL0330 

3VEL0340 

3VEL03b0 

3VEL0360 

3VEL0370 

3VEL03b0 

3VEL0390 

3VEL0400 

3VEL0410 

3VELn420 

3VEUD430 

3VEL0440 

3VEL0450 

3VEL0460 

3VEL0470 

3VEL0480 

3VEL0490 

3VEL0S0Q 

SVeLOSIO 

3VEL0520 

3VeL0Ij30 

3VEL0540 

3VEU0530 
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UUUOSo 
OOuO'j? 
UUUObii 
UUj0b9 
OUuOucI 
OUUOul 
UOUUo?. 
000063 
OUU06‘4 
00006b 
OOOOuo 
0 0 0 0 o 7 
000066 
00O069 
000070 
000071 
0oo07*d 
0UU073 
000074 
000U7b 
000076 
000077 
00O07& 
00O079 
OOOOoO 
000081 
OooOfic 
00O083 
000084 


XCK = U{I*1) 

YCK = V{IH) 

IK ,<XCK.LT.O.n) GO TO 110 

UCOMVA = -(4.»U(I + l)*(U(Itl)-U<I) )+U(l + lU(V(IIO-*V(IH+U'*V(l) 

i-V(I + l))+Sl»<U<I + l)-U(I”IJCPri*l) )/2,)*FACC 
GO TO lib 

UCOIJVA = -(4.*U(Ul)*(U(I+2)-U(Ul) )-»U(T + l)*{V(IH)+V(IH+n-Va) 
l-V(l + l) ) ♦Sl*(UtIH-fl)-U(I + l))/2. )*FACC 
COl-iTIIJUE 

IF (TCK.LT.O.O) GO TO 120 

Vk.OlJVA = -(U.»V<IH)*(V(III)-V(I))+V(IH)*{U(Ni).fU(IH+l>- 
lU (1 ) -U ( I H) ) +S2 ♦ ( V ( IH > -V ( I ri-i ) ) /2 . ) *FACD 
GO TO 125 

VCOWVA = -(4,»V{IH)*{V(IH+NCPR)-V{IH))+V{IH)*<U(IH+I)f 
lU(I + l)-Ua)-U(ItO )+S2*TV(IH+l)-V(It!) )/2.)*FACD 
COi-JTIiJUE 

Ir ( (lOUIT-lSTRT) .EO.O) UCONVA = 0. 

AHl(I) = UCOtJVA ’ 

At12(I) = VCONVA 
CCtwTIrJUE 

IF (lORIG.EQ.D) GO TO 200 
F = A.‘1llYm/(D(n»*0.3333l 

Gi = l.+F*OELT*( (lo.»U(l+l)»U{I + l)+Sl*SlJ»»0.b)/( (D ( I ) •* D ( H 1 ) ) **2 
A1 = U(I-tl)-tFAC*(D(I)+D(I+l) )*tH(n-H(I+T) )+X»DEl-T+FAA«Sl*COR 

1+UCONVA 

G2 = l.+F»yCLT* ( (l6.*V(IH)»V(IH)+52tS2)**0.5)/( (D(I)+D(IH) )»*2) 

A2 = V(I|i)+FAC*tD(I)+D(lHn*(HtI)-H(IH))+y»:DELT-FAA*S2»C0R-‘VC0NVA 
UP(I+1) = A1/|G1*FX( I+l) ) 

VP(IH) r a2/(G2*FY(IH) ) 


0000o5 

005 


GO TO 300 

00008b 

005 

200 

COtlTIWUE 

0U0U87 

005 


F = AMtlVm 

000088 

UOb 


UP( I+l) = U(I+1)+FAA4S1»G0R+FAC*(D(I)+0(I+1> )*(H(I) 

000089 

005 


l-H(I + l) )+X*DELT-F*DELT*U(l4l)’YABS(U(l+-l) )/< ( tOd) 

000090 

005 


2+01 1+1) 1/2. )*»2i3333) 

000091 

005 


VeMIll) = V(.IH)+FAA+S2*C0R+FAC+(D<I)+D(IHj ) ♦ (H ( I ) “H ( IN) ) 

000092 

1105 


1+ Y*DELT-F+dVlt«V( IH)+A nS( Y( IH) )/( ( (0( I)+0(IH) )/2. ) **2.3333) 

000093 

005 

300 

COUTINUE 

00u094 

U07 


IF ( (lOUlT-ISTf(T) .IJE.O) 60 TO 310 

000095 

005 


UP(IQUIT) =0. 

OOoOyo 

005 

305 

UP(1(5UIT+1) =0. 

000097 

005 

310 

CONTIMUE 

000090 

005 

cc 


000099 

005 

350 

LL = J 

000100 

005 

355 

ISTRT = IHI'IPL (LU 

OOOlOl 

005 


IGUIT = IBWDRU-1-) 

000102 

005 


DU 380 I = ISTRTrlOUIT 

0UO103 

OUb 


IF ((1STHT-IOUIT).EO.O) GO TO 360 

00O1O4 

■ 005 


F = AMnY(I)/(D(n**0.3333) 

OOOlOb 

005 


SI = V(I) + V(I+1) 

OOOlOb 

005 


Gl = 1.+F*0ELT* ( (16.*U(I+-1)*U{I+1)+S1*S1)»*0.5)/ 

000107 

005 


1 ( (5( I ) +D( I+l ) ) ‘*2) 

OOOIOU 

005 


Ai = U(I + l)+FAC*(3(I)+DU + l))*()H(I)-H(I+-l))+X*DEtT 

00O109 

005 


1+FAA*S1*C0K 

000110 

005 


U)-'( I + l) = Al/Gl 

uuoiii 

005 


GO TO 380 

000112 

005 

360 

UP(I+1) = 0. 


3VEL0560 

3VEU0570 

3VEl-05a0 

3VEL0b90 

3VEL0600 

3VEL0610 

3VEL0620 

3VEL0630 

3VEL0640 

3VeL0650 

3VEL0660 

3VEL0670 

3VEL0680 

3VEL0690 

3VEL0700 

3VEL0710 

3VEL0720 

3VEL0730 

3VEL0740 

3VEL0750 

3VEL076O 

3VEU0770 

13VEL07QO 

3VEL0790 

3VEL0800 

3VEL0810 

3VEL0820 

3VELOS30 

3VEL0840 

3VELU850 

3VEL0860 

3VEL.0870 

3VEL0860 

3VEL0890 

3VEL0900 

3VEL0910 

3VEL0920 

3VEL0930 

3VEL0940 

3VEL095C 

3VEL0960 

3\/EL0970 

3VELf)9aO 

3VEL0990 

3VEL10Q0 

3VEL1010 

3VEL1020 

3VEL1030 

3VEL1040 

3VEL1050 

3VEL1060 

3VEL1U70 

3VEL1060 

3VEL1Q90 

3VEL1100 

3VE1.1110 

3VEL1120 



90 


ociom OUi) WH(I(i01T) =0. 

oojim 00b jou coinioiif 

OUUUb OUb LL = LL+1 

UOUUO 007 IF (LL.LC.IbBl GO TO 355 

0UuU7 UUd C RiVEH !-‘AH5H D. 

OOolld OOd C 

000119 00b C 

OUU120 UOb U' {IFR0.v.tC.'.99999l GO TO “490 

000121 00a DO 450 K = 1 .2 

000122 005 I5T:a = IRf'li.-'nirO 

000123 005 luUIT = lRMrT(t', ) 

000124 005 DO 450 I = ISMiTr IOUlT,fJCPR 

U00125 005 IF ( ( I'tUCPK) .GT.tiC) 60 TO 390 

OOU12u 005 51 = V(I) + Vtl + i) + Vd+NCPR) + V ( I+NCPR-*-! ) 

000127 005 60 TO 395 

00U120 005 390 Si = Mil) t-V(I + l) ’ 

000129 005 395 mCKT:H<I) - H(l+i) 

000130 005 KtK2 r: 2r.il - H(I+U 

00U131 UOb HtK3 = H(I) - 2RB 

00U132 005 Ip (HCKl) 401,400,402 

000133 005 402 Ir (KCK3) 400,400,410 

000134 005 401 IF (HCK2) 410,400,400 

OU0135 005 410 DU = 0 . 5f (li ( I + l MH ( I ) ) -ZRB 

O0ol3o OUb Ai = U(I + l)4FACT(0n)+D(I + l) )*(H(I)-H(I+1) )4-X*DEU+FAfl»5l*C0R 

000137 007 61 = 1.4DEI.T» (L,(1 + 1)+D(I) ) «AUS(U(I + 1) )/(0.18*OB»DU*DELS) 

0U0138 005 UP(l+l) = Al/61 

000159 005 60 TO 450 

000140 005 400 UPtl+1) =0,0 

00014.1 005 450 CO(,TI.MU£ 

000142 005 490 L'O 50U I = ItIV 

00U143 005 J = lUC(I) 

OOOlH'f 005 500 VP(J + DCPH) = 0.0 

0Dul45 005 CC N'O HEhb TO 2ER0 THE VELOCITIES OF THE RIGHT MOST CELL 

000149 005 CC OH EACH SEG^EHT 

000147 005 DO 510 I = 1»IU 

000141 005 >0 = ImC'U) 

000149 005 510 UP(J+1) =0,0 

UUOL50 005 return 

00U151 005 EUD 

END ELT. 


3VEL1130 

3VEL1140 

3VEL1150 

3VEH160 

3VEL1170 

3VELU80 

3VELU90 

3VEL12D0 

3VEL1210 

3VEL1220 

3VEL1230 

3VEL1240 

3VEL1250 

3VEL1260 

3VEL1270 

3VEL1280 

3VEL1290 

3VEL1300 

3VEL1310 

3VEL1320 

3VEL1330 

3VEL1340 

3VtL.1350 

3VEL1360 

3VEH370 

3VEL1380 

3VEL1390 

3VEL140D 

3VEL1410 

3VEL1420 

3VEL1430 

3VEL1440 

3VEL1450 

3VEL1460 

3VEL1470 

3VEL1480 

3VEL1490 

3VEL1500 

3VEL151Q 
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iillLUTtL SAi^«ijUbH 

tHUT7 KLib70 ul/*9-l7:46;i5-(28f ) 


000001 

02b 

C 

00OOU2 

U2tt 

C 

uouOOi. 

U2b 

e 

UOOU04 

02u 

c 

OOoUOb 

02d 

c 

UOuOOo 

024 


00O0U7 

024 


uUuuOd 

024 


OOUU09 

024 


OOuOlO 

024 


UOOOll 

024 


OUUOlB 

02t 


0U0013 

024 


UOOU14 • 

024 


0UuO15 

024 

10 

OOoUio 

024 


000017 

024 

11 

0 0 U ti 1 b 

024 


ooooiy 

024 

20 

00UO20 

024 


0U0021 

024 

21 

00UO22 

024 


0UOU23 

024 


000024 

024 


0UJO2a 

024 


00002O 

024 


00(Jt)27 

024 

ccc 

OOU02B 

024 

cc 

Q0U029 

024 

c 

OUUO3O 

024 

3G 

0UDU31 

024 


0OU052 

024 


OUU033 

024 


OO0U34 

024 


000033 

024 


000036 

024 


000037 

024 


000036 

024 

50 

U0U63'4 

024 


00UU4U 

024 


O00041 

024 


OUU042 

024 


000043 

024 


000044 

024 


00oO4‘j 

024 


00u04u 

024 


000047 

024 

70 

U0U04d 

024 


000049 

024 

BO 

OOuObO 

024 

ioo 

UOU051 

024 


000052 

024 



VtKSIOr; 1.2ft. 7c 

♦ 

SUiJKOUTlNE 5BA!!Kl(IllP.lP*NCPR.DELT,r.ELStFAA.PACfC0RfX»Y) 
OiMENblOiJ IMS5C(100)» 2CHU00). lAXlS(inO) 

COlWOil /A/ 7(2£iU0). H(2500). D(2500). lBnUU(2UCJ), IUNDRC200)# 
1NU,MIU0). 1REP(200) 

COmmO’J /»')/ VMACv(2500), BETA (2500). AHU2500). AM2(2500)r 
1SJR|X(2500) . SUI-iY<2500) . SOMXNT { 2500 ) . SUMYnT (2500) , 

20(2500). V(2500). UP(2500). VP(2500) 

INTEGER. 1EST(2) /'X'.'T*/ 

R£A0( HIP. 10) M5dC 
FORMAT (1 5) 

VRlTEdP.m MSBC 

format (• '.151 ^ 

REAO( lUP. ?,0 ) ( INSbC ( 1) . TAX IS ( I ) .ZCH( I) . I = 1.I'(5BC) 
FjRMATt5(Ii5'Al.F'+.l> ) 

(\KITE( IP.21) (IlJS.iCd) »1AXIS(I) .ZCH(I) .I=1.NSBC) 

FOKmaTC • (Bt J5rAl.F4.in 

EiJTRY SBANK2( IIjP . IP.NCPR r DELT.DEUS.FAA.FAC.COR. X. Y ) 

Du 100 J = I.NSbC 
I =■ lUSOC(J) 

IF (lAXIS.EO.TESTd)) 60 TO 50 
IF (IaXIS,E0.TEST(2)) GO TO 30 

CALL ERROR! ) 

D3 = 0.5*(H(I)-»mm))-ZCH(I) 

51 = V(l)i\/(I+1)+V(MCPR+I)+V(NCPR+I+1) 

IFJDB.LT.O.O) GO TO 70 

61 = 1,+DELT* (D( I + II+DCI) ) ♦AB5(U(I + 1) )/'!C. ld*DU*DB»DELS) 

Al = U(I + i)+FAC»(D(n+0(I + in*tH(I)-H<I + l) ) 

1+A»DELT+FAA*SI«C0r 
UP(I+ 1) = Al/61 
GO TO 100 

C)I3 = 0.5+(1i(I)+H(I+ITCPK))-ZCH(1) 

52 = u(l)+U!l + n+U(nCPR+I)+U(NCPR+l+l) 

IF (DB.LT. 0.0) GO TO BO 

g 2 = i.4DELT*(D(I+i'lCPR)+D(I) )*AB5{V(I+WCPR))/ 
1(0.1B»DB»DB»DELS) 

A2 = V(I+NCPR)+FAC+(D(1)+D(I+NCPR) )»(H(I)-H(ItNCPR) ) 
l + Y*DLLTt-FAA.S2*C0R 

VP(I+NCPR) =. a 2/G2 
GO TO 100 
UP(Itl) = 0.0 
GO TO 100 

VP(I+|JCPR) = 0.0 , 

COMTIilUE 
RETURN 
, E.O 


ASrKOOIO 

uS[)K0020 

4SbK0030 

4SBK0040 

4SBK0050 

ASnKOOfaO 

HSBK0070 

4SUKOOBO 

4SBK0090 

4SHKOIOO 

BSbKOIIO 

45rK0120 

4SdK0130 

4SBK0140 

45BK0150 

4SUKU160 

4SBK0170 

45BK01B0 

4SOK0190 

4SBKC200 

4SBK021O 

4 S 13 KO 22 O 

4SBK0230 

m5bK0240 

4SBK0250 

45BK02fa0 

4SBK0270 

4SHK02d0 

4SBK0290 

4SBK0300 

4SBKO3IO 

4SBK0320 

4SBK0330 

4SBK0340 

45(3X0350 

4SBK0360 

4SbKQ370 

4SRK03HO 

4SBK0390 

4SbKOhOO 

4SBK0410 

4SBK0420 

4SRK0430 

45bK04MO 

4SBXO4S0 

4SBK0460 

45RK0470 

4SrK04B0 

45(3X0490 

4SDK05DO 

4SbKu5I0 

45pK0520 


END ELT. 
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ulILI rL LiAi'i.^iULia 
fcLIUT? uLiii7u 'Ji/ 


OUOUUl 

UOuCiUi:' 

OOUOOJ 

OOuOUl* 

(JUUUOb 

UUuUOb 

00U007 

OOUUOo 

0UU009 

OOUOlO 

OOoOii 

oouuia 

00UUi3 

oodum 

uuudlj 

uuudlo 

0 Odd 17 

00u0i6 

00dU19 

Uguoau 

U0UO21 

udooaa 

oouuao 

U00024 

00UU2U 

00o02t> 

0UO0E7 

0UU023 

U0dU29 

OOU03O 

UUdOJi. 

OJU032 

OOU033 

UOu03‘+ 

00d055 

lJOOU3o 

00d037 

OdOdOd 

00U039 

0000*40 

OOdOol 

000042 

000043 

000044 

00d04b 

0U004ffi 

OOUJiK' 


024 

024 

024 

024 

024 

020 

020 

020 

020 

02d 

020 

020 

020 

020 

020 

02d 

d20 

020 

020 

020 

020 

020 

020 

U20 

020 

02d 

020 

020 

020 

020 

d2d 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

020 

U20 

020 


29-l7:4£»;ig-(24f ) 

c 

C 

c 
c 
c 


cc 

cc 


Vu^OIOI. 1.2Q.7ti 
♦ 


SOimOuTIilE MA6U'lR(I5BMJCr>niNC,AHRAD.Pn 
qlReCTIOIJ 15 MCaSURED from X-aXIS COUMTER-CLOCKWISE 

COMMON /A/ 2(25)00), H(2500)r 0(2500) f lBNDLt200), IliNDK(200), 

^CUt-’plOiT^/U/^VM^ BETA(2500). AHK2500), AH2(2500), 

15UM(2500), SUMy(2b(IO), SUMXNT ( 2500 ) . SUMYNT ( 2500 ) , 

2U(2500), V(2500)* UP(2500), VP(2500) 

DO 150 K = l,.I5ri 

Ir (ItiliOUtK) .GT, (NC-UCPR) ) GO TO 160 
13TRT = ItjNDUK) 
luUlT = InhrR(K) 

DO lOd I = ISTRTrlUUtT 
IF { (I0Un-15TRT),E0,0) 

Go TO 40 


30 


40 

50 

ud 

70 


£Jd 


CC 

at 

do 


av 


90 


GO TO 30 


VV = V(I) 

Ud = 0. 

GO TO 50 

Ud = tU(I)+U(I + U)/2. 

VV = (V(l3+V( I+NCPR) )/2, 

IF (UU) 70#ft5,e0 
IF (VV) aO,9D»90 

BtTA(I) = {ATAN(W/UU)+PI)vANRAD 
Vi'lAG(I ) = ( (UIJ*UU+VV*VV) ♦♦0 ,5) /D( I ) 
GO TO 100 

UtiTAU) = (ATAN(VV/UU)+2.*PI)»ANRAD 
V,-iAG(I) = ( (UU»UU+VV*VV)* + 0.5)/D(I) 
GO TO 100 
CHECK above 
IF(VV)B 6, 100,87 
BcTA(l) = 270, 

Vi'i«G(I) = AB5(VV)/D(1) 

60 TO 100 

beta ( n = 

VKAG(I) = 

00 TO 100 
pETAin = 

VMAG(I) = 

100 CONTI, A)E 

iBO Continue 
lt,0 KETURil 
E«U 


90. 

VV/D(I) 

<ATAN(VV/UU) )»ANRAp 
( CUU*UUtVV+VV) »*Q.5)/D(I) 


bHGOOOlQ 

bMGD0020 

5MGUO030 

bMGD0040 

bMGDQOBO 

bMGDOOoO 

bMGCi0070 

BMGDOOao 

bMGD009Q 

bMGDOlOO 

5MGU0110 

bMGDOiaO 

5MGDU130 

bMGn0140 

5MGD0150 

bMGD0160 

SM6DU170 

bMGLidiaO 

5M6CD190 

5MGD0200 

5MGD0210 

5MGD0220 

SMGD0230 

5MGP0240 

5MGC0250 

5MGD0260 

5M6D027Q 

bMGD02B0 

5MGD0290 

5M6D0300 

bMGD0310 

bMG00320 

5MGD0330 

5MGD0340 

5MGD03bO 

bMGD0360 

bMG00370 

5MG00380 

5^600390 

5MGD0400 

bM6004iO 

5MGC042O 

5MGDO430 

5MGD0440 

bMGBa450 

5MG00460 

5MGD0470 


END ELT 
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ijel r f l 

bAi'i.SUUb 



ECT0T7 

HLib70 01/29-17 

;4f);24-(59. ) 

UOUUOl 

059 

C 


U0uU02 

059 

c 

♦ 

OUU003 

059 

c 

VERSION 1.28.76 

0UU004 

059 

c 

4 - 

UUO0J5 

059 

c 

SU5R0JtliJE DISCMldl.'P, IP, DELS. DELT , VOL , DI5R ,SD ISRp, SOIbP.tJ.LUSCH) 

0OUOO6 

055 


00OU07 

055 


L'Ii/iLM510II uISR( 20) ,5DISRP(20) ,SDlSRfJ{20) ,VOL(20) ,UCDSCH(20) r 

000008 

Ob5 


lICl}SCil(2U,5i)), IAXIS(2ni 

ooouoy 

055 


COMlOil /A/ Z(2500), H(250U). D(2500). IdNDL(200). IUNDR(2uO). 

OOoUlO 

055 


lt,’J'<H100),IREPt200) 

oouou 

055 


CbM..10>J /IV VMAb(2500). »ETA(2500). AH1I2500), AH2{2500), 

OOU012 

055 


lSuWX(25(rO) . 5ur-iY(2j00) » SUMXIJT ( 2500 ) . SUMYNT < 2500 ) , 

0UUU13 

055 


2U(2500>, V(2500), UP(2500). VP(2bOO) 

OUU014 

055 


INTE6LI! IT£ST(2) /'X'.'Y’/ 

UOviDlb 

055 


reach IliP. 10) LD5CH 

UUWOlo 

055 

10 

FOriNATllb) 

000017 

05b 


WiUTEdP.lU LDSCH 

OOuOlo 

055 

11 

formats *.15) 

000019 

055 


DO 15 I = 1, LDSCH 

000020 

055 


DISR<1) =0.0 

00u021 

055 


VuL(l) = 0.0 

OUU022 

055 


SDISRPII) = 0.0 

000023 

055 

15 

SDISRim =0.0 

000024 

055 


REAOaHP.20) (NCDSCH(K). 1AXIS<K), K=1 .LdSCH) 

UOU025 

055 

20 

FbRMATllb( 14, Al ) ) 

000026 

055 


lv2iTE(lP.?l) (HCDSCH(K), lAXIS(K), K=1 ,LdSCH) 

00U027 

055 

21 

format (• «, l6tl4,Al)) 

000028 

05b 


DO -2C J = l.LDbCH 

000029 

055 


JJ = NCOSCH(J) 

000U30 

055 


READ ( U!P, 25) (1CD5CH(K,J),K=1.JJ) 

000031 

055 

25 

FuR 'iAT (1615) 

0U0032 

055 

26 

KfaTE(lP,2o) (ICDSCtHK.Jl 

000033 

055 

26 

FORMAT (< <,16I5) 

UUU034 

055 


ENTRY D15CH2(lliPrIP,DELSrUELT.VOL,OlSR,SDlSRP,SDISRN,L05CH) 

OOU035 

055 


DO lOu I = 1, LDSCH 

OPJi'036 

055 


TbTAUV =0, 

000037 

05b 


TOTALu =0. 

000038 

055 


NK = NCDSCrUI) 

00UU39 

05b 


If (lAXlScI) .EL.ITESTU) ) GO TO 30 

000040 

055 


IF tlAXISd) .E6.TTE5T(2)) GO TO 40 

000041 

055 

C 

Call error ( ) 

000042 

055 

30 

DO 50 J = l.NK 

00U043 

055 


Kk = ICDSCHCJ,!) 

000044 

055 


TOTALu = U(KK) + TOTaLU + UP(KK‘- 

0U0O45 

055 

50 

CONTTNHE 

00004b 

058 


DiSRd) = DISR(I) + TOTaLU/2. 

000047 

055 


GO TO 100 

004048 

055 

40 

DO oO J = 1,NK 

0U0U49 

055 


KK = ICDSCIH J, I ) 

UOuUbO 

055 


TOTaLV = V(KK) + TOTaLV + VP(KK) 

UOUUSl 

055 

oO 

CuNTTiaUE 

0UOO52 

055 


DlSRdi = DISRd) + TOTALV/2. 

0QUU53 

055 

lUO 

CbNTIwUE 

000054 

055 


return , . 

000055 

055 


ENTRY 0I5CH3( INPfllH0EL5,DELT,V0L,DISR.5D15RP,SDlSRN, LDSCH) 

CiOuObo 

055 


DO 150 I = IfLbSCH 

UUUiJb ? 

U5d 


VOLd):= Dl5Pd)*DELS*T)ELT/3600, 

0OJ058 

055 


Ik (VOLd)) 110,111,112 

OOU059 

Ci55 

11 0 

bblSRMdl = SDlSRlHI) + VOLdl436ao. 

UU0060 

055 


(,U TO 111 

OuoUoi 

050 

112 

SDlSRHd) = SDlSRPdl + VOUIl+3600. 

00OU62 

055 

111 

CUNT UlUE 

000063 

055 


DlSHtn - o.t' 

00U06-4 

055 

150 

CONTINUE 

OOoOoa 

055 


return 

OOoUbb 

055 


END: 


gdscooio 

6D5C0020 

(,DSC0030 

bDSCOOUO 

6DSC00‘j0 

6DSC006O 

6DSC0070 

6DSCOU80 

6DSC0090 

bDSCUlOO 

6DSC0110 

6DSC0120 

6DSC0130 

6DSC0140 

6DSC015O 

6D5C0160 

6DSC0170 

6DSCU180 

60SC0190 

60SC0200 

60SC0210 

6D5C0220 

6DSC023O 

60500240 

6DSC0250 

6DSCQ260 

6DSCC)270 

60500260 

60SC0290 

60S00300 

6DSC0310 

60500320 

60500330 

60SCD340 

60500350 

60500360 

6DSOU370 

60500360 

60500390 

60500400 

6DSC0410 

60500420 

60500430 

60500440 

6050 0^ 50 

60500460 

60500470 

bOSCOHBO 

60500490 

60500500 

6DS00510 

60500520 

6DSC0530 

60500540 

60500550 

6DSC0560 

60500570 

6050 0580 

60500590 

60500600 

60500610 

60500620 

6DSCU630 

60500640 

6DSC0650 

60500660 


ENu ELT 
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V 


oouuui 

046 

C 

UUUU02 

04o 

c 

i)0o0(J3 

U4d 

c 

0000U4 

04u 

c 

000005 

U4(> 

c 

OOOOOo 

042 


00U007 

042 


00(;0yd 

042 


ooojjy 

042 


(lOuOlO 

042 


OOUUll 

042 


00UU12 

042 


0000i3 

042 


000014 

042 


OOJUlb 

042 


OOUUIu 

042 

c 

000017 

042 

c 

OUUOld 

042 

c 

000019 

042 

ccc 

OOOU20 

042 

cc 

000021 

042 

c 

000022 

042 


0UuC23 

042 

100 

000024 

042 


000025 

042 

101 

00 00 2d 

042 


000027 

042 

c 

00o023 

042 


000029 

042 


000030 

042 

120 

000031 

042 


000032 

042 

121 

000033 

042 


000034 

042 

150 

0uo035 

042 

c 

OOOOiu 

042 


000037 

042 

160 

000036 

042 


0000,59 

042 

lei 

OOOU40 

042 

c 

000041 

042 


000042 

042 

c 

000043 

042 


000044 

042 

300 

O0O045. 

042 


000049 

042 

301 

000047 

042 

C 

000040 

042 

c 

000049 

042 


00O050 

042 


000051 

042 

3u5 

000052 

042 


000053 

042 

3u6 

000054 

042 


OU0055 

042 
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VuHSION 1.26.76 
♦ 

5U3f!0UTII'jE ROUMDl ( INP. IP.TMIL. ITSMiTD.HDB.FAB.FaC.DELTiDELS. 
lUCPR.fO 

DlWEHSlOH HDB(20) rTD(20),CA(20».CD<20) -CC(20) rCD(20». 
lHuC(2D) .(1ATURC(20> ,UlliF(2U) .II'4RCt20*50) .ZDt20»60) 

CJi'iMO’! /A/ Z(2b0n), H(2SrjO). D <2500 ) r lBNOL( 200 ) . l3NDH1200)t 
iMUMdOO). IPEP(200) 

COf^MO.-J /B/ VMAO(2500)» 3ETA(2500>. AH1(2500). AH2(2bOO). 

16JMX t250U). SUf'Y(2500). SuMXtJT ( 2600) . SUMYNT (2500 ) , 

2U(2500). V(250U)» UP{2500)* 0P(2500) _ . 

lUTEGtlP TEST(6 ) / • 5U» i » S3' » »SR * . ’ 5T ’ . 'RL* . ’ RE ' • ’ RR* • »RT»/ 
fcjOUlMDARY CELLS ARE THOSE CELLS NEST TO THE LEFT , RIGHT. 

TUP, uR BOTTOM MOST CELLS TM EACH ROW OR COLOUMN 


« OF tide station will be used 


read ili boundary conditions 
initial calculations 
READ III starting TIME & 

READ ( IliP . 100 ) TMIL. ITSN 
format (F b.2. 15) 

WxlTEilP.lOl) TMIL.ITSN 
FOHi.1AT(t f.F5.2.l5) 

READ il/cOEF5°AND TIME PHASE USED TO CAL THE TIDE FUNCTIONS 
DO 150 I = 1. ITSM 

READ(INP,120) CA( I) .CB(I) .CC(I) ,CD(1) »TD{1) 

FOHMAT(5F10,6) , 

WKITEdP. 121) cat I) .CBd) .CC(I) »CD(1) »TD(I) 

FOnTAId ’ fS-FlO.G) 

TUd) = TDd) + 0.5*DELT/60. 

HDRd ) = CA ( 1 )+CB( 1 ) ♦COS (CC d ) »TD( I ) +CDd) ) 

read total H of boundaries not including RIVER MARSH BOUNDARIES 

READdNPdoO) NTB 

FORMATdS) 

wi-tlTEdP.iai) ntb 

FOkMATI* '.15) 

SPECIFICATION OF BOUNDARIES 
DO 400 I = I.NIB 

RtAD IN fl OF CELLS IN THE BOUNDARY. ITS NATURE. AND INF 
REAUdNP.300) NATURE(I) .NBCd) .RINF(I) 

FORMAT( A2. I3.F1U.3) 

wRITE( 1P.301) NATURE(I) .NbCCI) .RiNFTi) 

FORMATC '.A2.T3.F10.3) 

EXPLAIN ABOVE 
read in CELL «S 
NK = NRCd ) 

READdNP.30D) ( INoCd.U) fU=I.NK) 

FUKMATdblS) 

Wt<irEdP.30C) (IHUCd.J)^U=l.NK) 

format (•» .1615 > 

IF ( (naTUrEI I ) .EG.TCST dl ) , OR. (NATURE t I ) .EQ,TEST(2 ) ) ) GO TO 307 
GO TO 331 


7BND0010 

7BND0020 

7BND0030 

7BND0040 

7BND0050 

78ND0060 

7BND0070 

7BND0060 

7BND0090 

7BND0100 

7BND0110 

7BND0120 

7DND0130 

7DND0140 

7BMD0160 

7BND0160 

7DND0170 

7BND0160 

7BND0190 

7BND020O 

7BND0210 

7BND0220 

7BND0230 

7BND0240 

7BND0250 

7BND0260 

7DND0270 

7BMD0280 

7BND0290 

7BND0300 

7BND0310 

7BNDC320 

7BND0330 

7BND0340 

7BND0350 

7BND0360 

7DND0370 

7BND0380 

7BND0390 

7BND0400 

7BND0410 

7BND0420 

7BN00430 

7BND0440 

7bN00450 

7BND0460 

7BN00470 

7BND0480 

7BND0490 

7BND0500 

7BNDU510 

7BND0520 

7DND0530 

7BMD0540 

7BND055Q 
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OOUObu 

U<42 

3j7 

^(i£:;*,r)uii(’»30B,) (Zbii»k)* k = IiNK) 

00U0B7 

042 

308 

roH'.iArtioF5.o) 

OOOCIBB 

UH2 


WHITE ( IB t 305) (ZDtliK)* K = 1»NK) 

00UUB9 

042 

309 

FOI'.MMC* 'fibFB.O) 

UU JU6U 

042 

310 

IF (riATUUEID.LO.TESTd)) GO TO 315 

U0UU61 

042 


GO TO 321 

OOu Ju2 

042 

315 

DO 320 K = 1»NK 

0UU063 

042 


WA = ILHCtl.K) 

00UU64 

042 


JJ = (UNK(l) 

OuUOoB 

042 


UH(t<K) = U(KK)+FAC*(D(KK)+HDB(JJ)+ZB<IfK) )*(hDB(JJ)-H(KK) ) 

OUUObu 

042 

320 

COliTIfjUE 

000 U67 

042 


GO TO 400 

OOOUub 

042 

321 

IF (NATUflEin.CG.TESTlZ)) 60 TO 325 

U0UU69 

042 


GO TO 331 

00U070 

042 

325 

Ou 33vi K = 1»NK 

000071 

042 


KK. r 1 i;8C(I»K) 

C0UU72 

04e. 


JJ = KINF(I) 

000073 

042 

330 

VP(KK) = VIKKMFAC»(DIKK)+HDBI JJ)+ZBtI»K) )*(HDBlOJ)-H!KKn 

00UG74 

042 


GO TO 400 

000075 

04^ 

331 

IF (NaTUKE(I).LQ.TL'ST(3)) go to 335 

000070 

042 


GO TO 341 

000077 

042 

C 

RIGHT MOST ROUHDARY 

000070 

042 

335 

DO 340 K = l.NK 

00U079 

042 


KK = IMBCUrK) 

000080 

042 


JJ = RINFM) 

000081 

042 

340 

H(KK) = hDB(JJ) 

000082 

042 


60 TO 400 

00U0O3 

042 

341 

IF (NATUr<Em.E0.TEST(4)) GOTO 345 

0UU064 

042 


GO TO 351 

000085 

042 

C 

TOP boundary 

0UUU8O 

042 

345 

DO 350 K = 1*NK 

000087 

042 


KK = IhRCa.K) 

OOOOOo - 

04t 


J^J z R I^jF ( I ) 

0uoUu9 

uh2 

350 

hlKK) = HDB(JJ) 

000090 

042 


GO TO 400 

00O091 

042 

351 

IF UvaTUKc(I) .La.TEST(b) ) GO TO 355 

000092 

042 


GU TO 361 

000093 

042 

C 

RivcK Entrance from left 

0OO09H 

U4<; 

055 

DO 360 K = 1»NK 

000095 

042 


KK = INBCtl.K) 

00009b 

042 


UP(KK) = RlNF(I) 

U0u097 

042 

3u0 

UlKK) = RINF(I) 

000098 

042 


GO TO 400 

000099 

042 

3b 1 

IF (|JaTURE(I).EG.T£ST(6)) GO TO 365 

000100 

042 


GO TO 371 

OOOiUl 

042 

3u5 

Dv/ 570 K = liNK 

000102 

042 


Kiv = innClIfK) 

000103 

042 


V('(AK) = HINFU) 

000104 

042 

370 

V(KK> = rUMFtI) 

OOUiOb 

042 


GO TO 400 

OOOlOb 

042 

371 

IF(UATURE(X) ,EQ.TESTt7)> GO TO 375 

000107 

042 


C-o TO 381 

000108 

042 

375 

DO Jeo K = 1 »NK 

000109 

042 


KK = ItlHC(I.K) 

OOullO 

044 


f!j = -RINF(I) 

OOulll 

042 

380 

HlKK) = |itKK)+FAB*(UPlKK)+VP(KK)-KB-VP(*'-K4UCPR))+H*UELT 

000112 

042 


60 TO 400 


7BlJDOb60 

7BND0b70 

7BND05B0 

7BtiD0590 

7BND0600 

7BWD0610 

7BMD0620 

7BND0630 

7BMD0b40 

7BND06B0 

7BND0660 

7BrjD0670 

7BUD0680 

713^00690 

7BND0700 

7BND0710 

7BND(j720 

7BI1D0730 

7BtJD0740 

7BND07B0 

7BND0760 

7BU00770 

7BND0780 

76tJD0790 

7BMD0800 

7BND0610 

7BND0820 

7BND0630 

7BND0840 

7BND0a50 

7BND0860 

7nMD0870 

7BND0880 

7BND0890 

7BWD0900 

7BNB0910 

7BND0920 

7BND0930 

7BUD0940 

7BND0950 

7DND0960 

7BND0970 

7BND0960 

7BNCJ0990 

7BfJD1000 

7BN01010 

7DND1020 

7BHD1030 

7BND1.040 

7BND1050 

7DND1060 

7BND1070 

7BND1080 

7BND1090 

7ONDU00 

7BND1110 

7BND1120 


reproducibility op the 
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0UU113 

04c 

3tU 

IF (liATUFEtn.BO.TESKan GO TO 385 

JUul.14 

042 

C. 

CmLL BI'P.OIU ) 

OUMUa 

042 

3a5 

[lO 39<l K = l.NK 

OUUUo 

042 


KK ,= II,-4CU»K) 

fJOUll 7 

044 


hJ z -i<I!iF(I) 

OOuilB 

042 

390 

H(KK) = H(KK)4FaB»(UP(KK)4VP(KK)-UP{KK+1)-RU)+R*0ELT 

0uoil9 

U4c 

4dO 

co.jTiijue 

UOiUBO 

042 


RETIJRU 

0UU121 

042 


rUTHY 60iJNd2(II'.Pt:P»TMlL,IT5NrTr)',HDB»FAB*FACfDELTfOELS,NCPR. 

00U128 

042 

CC 

CALCULATE TID hElGiUS FOR T AND BELT 

000X23 

042 


DO 550 I = 1»1TSN 

00U124 

042 


TU<I> = TD(I) + DELT/60, 

00Ui2b 

042 

550 

HUij(I) = CA(I)+CUtI)*CCS<CC(I)»TD(I)4CDa)) 

000120 

, 042 . 


DO 700 I = IrNTB 

000127 

042 


UK = NBC ( I ) 

00012a 

042 


IF (NaTURE(I).EO.TEST(I)) go TO 615 

000129 

042 


GO TO 621 

000130 

042 

bl5 

IF UmK.GT.I) go To 616 

000131 

042 


KK s It.BCd.ll 

000132 

042 


Jj = RIMFa) 

000133 

042 


UIT<KK) = U(KK) + FAC*(D(KK> + HDP(JO) ZB ( I r 1 ) ) ♦ ( HOB (oJ) -H( 

00O134 

042 


GO TO 700 

00ul3b 

042 

fal8 

DO 620 K = 1»NK 

U0U13O 

042 


KK r ItiBCUiK) 

000137, 

042 


JJ = RINF(I) 

00013a 

042 


G1 = l.+DELT*(D(KK)+2B(i*K)+HDP,(JJ) )*ABS(U(KK) )/(18.*DEL5) 

OOU139 

042 


Al = U(KK)+FAC*(D(KK)+H0BIJJ)+ZB(I»K))*<HDB(J0)-H(KK)) 

0UO140 

042 

b20 

UP(KK) = Al/Gl 

00O141 

042 


Go TO 700 

000142 

042 

bcl 

IF ( nature ( n .£0. TEST (2) ) GO TO 625 

000143 

042 


GO TO 631 

000144 

042 

625 

DO 630 K = l.NK 

00 01 4b 

042 


KK = INBCtl.K) 

00014O 

042 


JJ r R1NF(1> 

000147 

042 

b30 

VP(KK) = V(KK)+FAC*(D(KK)+Hn8< JJ)+ZD(I.K> JJ)-H(KKn 

00014B 

042 


GO TO 700 .. 

000149 

042 

o31 

IF (NaTURE(I).lO,T£ST(3) ) 60 TO 635 

000 160 

042 


GU TO 641 

oooibi 

042 

C 

RIGHT MOST 

004162 

042 

o35 

DO 640 K = IflaK 

000163 

042 


Ka = INBC(I.K)+1 

00ul64 

042 


KJ = iNfiC(i.K) 

00015‘J 

042 


JJ = tUNF(l) 

OOOlSu 

042 


hlKK) = HDfKJJ) 

000157 

042 


G1 = i.+DELT4(C(Kj)4Z{KK)+H(KK) )*ABS(U(KK) )/(1«,»DELS) 

000163 

042 


Al = U{KK)4FAC4(D(KJ)+Z(KK)+HtKK))*(HtKJ)-HUB<JU)) 

000169 

042 

ohO 

UP(KK) = Al/Gl 

0U0160 

042 


GO TO 700 

U00161 

042 

U4l 

IF (NATURE(I) ,£Q,TEST(4) ) GO TO 645 

000162 

042 


GO TO 651 

000163 

042 

0 

TOP MuST 

000164 

042 

b45 

DO b5u K = 1. nK 

OOOlob 

04c 


KK = lIlHCd .K) 

00016a 

042 


JJ = lUNFd) 

OOOlrj? 

04 c 

650 

MUK) = HDi3( JJ) 

000166 

042 


GO TO 700 

000169 

04c 

651 

IF (NATUKEdl.EO.TESTCS)) 60 TO 655 


7BNLJU30 

7DND1140 

7BNDU50 

7BNCai60 

7BHDU70 

7BND1160 

7BND1190 

7BND1200 

7BND1210 

7BND1220 

7UfiD12J0 

7BND12‘40 

7BND1250 

7BND1260 

7BNB1270 

713IJD12B0 

7BUD1290 

7Br)D1300 

7BIJC;i3lO 

7HN31320 

7BND1330 

7BND1340 

7BtJD1350 

78NC)1360 

7BND1370 

7Bt)D1380 

7BND1390 

7BHD1400 

7BND1410 

7BND1420 

7QUD1430 

7BND1440 

7BND1450 

7BND1460 

7BND1470 

7BND1480 

7BND1490 

7BND1600 

7Bn015l0 

7BND1520 

7BND1530 

7BND1540 

7BIJ01550 

7BND1560 

7BND1S70 

7BHD15H0 

78ND1590 

7BND 1.600 

7BN01610 

7BMD1620 

78ND1630 

7BMD1640 

70ND1650 

7BND1660 

7BHDlb70 

7BND1660 

7BND1690 
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0UU170 

0h2 


DO TO i:(,\ 

00U171 

042 

ubb 

no 660 K = 1»NK 

0UU172 

042 


KK = ll.BCdiK) . 

00U173 

042 


UK(KK) = I'iXI'.Fd) 

UOUlYH 

042 

(j&O 

U(KK) r HIuFd) 

U0uX7B 

042 


GO TO 700 

0uul7b 

042 

iicX 

IF (NATURE(I).E0.TEST<6)) GO TO 665 

000X77 

042 


60 TO 67X 

000X78 

042 

065 

CO 670 K = 1(MK 

0OUX79 

042 


KK = lUBCdfK) 

OOuXSO 

042 


VI' (KK) = RiKFd) 

OOoXBX 

042 

o70 

V(KK) s RIiJFCI) 

000X82 

042 


GO TO 700 

OOOX(i3 

092 

o7l 

IF (NATUHEd ) .E0.TEST(7) ) GO TO 675 

OOOXM 

042 


GO TO 6ttX 

OOuXOl) 

042 

o75 

DO 680 K = IfNK 

OOoXUb 

042 


KK s INSCd.K) 

OOUX87 

044 


RB = -RliJFd) 

000X83 

042 

pOO 

H(KK) - H(KK)+r-A[3»(UP(KK)+VP(KK)-RB-VP<KK+HCPR) )+R*DEDT 

000X89 

042 


GO TO 700 

OOuX90 

042 

60X 

IF (NATUREd) .tG.TEST(Sn GO TO 685 

00OX9X 

042 

C 

CALL ERROR < > 

00OX92 

042 

085 

no 690 K = X»MK 

O0oX93 

042 


KK = ILBCd.K) 

000X94 

044 


R8 r -RiNFd) . 

000X95 

042 


H(KK) = K(KK)+FA0*(UP(KK)+VP(KK)-UP(KK+X)-RB)+R*DELT 

000X90 

042 

o90 

CONTIkUE 

000X97 

042 

700 

CONTINUE . 

000X98 

042 


RtrURN 

000X99 

042 


END 


END ELT. 


70MO17OO 

7HND1710 

7BND17aO 

7BND1730 

7BND:7‘*0 

7QDD1750 

7BND1760 

7BND1770 

7BND1700 

7BUD1790 

7BHD1600 

7BND1610 

7BIJD1820 

7BND1830 

7BND1B40 

7BND1650 

7BND1660 

7BfJD1670 

7Bt4D1880 

7BrJDi890 

7BNDI90Q 

7BND1910 

7BMD1920 

7BNDI930 

tBNDig^O 

7BtgU1950 

7BIJD1960 

7BND1970 

7BUD1930 

7BND1990 


OF 


f IS 500B 




